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fluorescence  technique  was  developed  for  time-resolved  measurements  of  species 
concentration  profiles.  This  technique  was  combined  with  laser  velocimetry  to  obtain 
measaieaamts  of  turbulent  mass  transport  in  a  flow  with  polymer  injection  and,  for 
comparison  purposes,  a  flow  with  water  injection.  Two-component  laser  velocimeter 
measurements  were  made  in  both  flows  to  examine  the  effect  of  both  the  injection 
process  and  the  evolving  polymer  concentration  field  on  the  structure  of  turbulence. 

In  the  flows  with  injection,  the  major  transport  mechanism  was  identified  as  long  (in 
the  streamwiae  direction)  filaments  of  high  concentration  fluid  periodically  lifting 
away  from  the  near-wall  region.  When  polymer  was  injected,  mixing  of  the  injected 
fluid  with  the  outer  flow  required  longer  streamwiae  distance  than  for  water 
injection.  The  high  extensional  viscosity  of  the  polymer  solution  resulted  in  a  strong 
correlation  between  large  instantaneous  polymer  concentrations  and  small  wall- 
normal  velocity  fluctuations  causing  lower  levels  of  turbulent  mass  transport.  After 
the  polymer  solution  began  to  mix  with  the  flow,  the  turbulent  mass  transput,  was 
smaller  than  in  the  flow  with  water  injection  but  the  degree  of  correlation  between 
velocity  and  concentration  fluctuations  was  comparable.  This  indicates  that  ttye  mass 
transport  processes  were  similar  for  these  two  cases,  hi  both  cases,  the  high  Schmidt 
number  resulted  in  high  concentration  fluid  mixing  from  the  near-wall  region  to  the 
outer  flow  and  then  back  near  the  wall  yielding  low  levels  of  net  turbulent  maw 
transport.  7 

Measurements  of  turbulence  structure  show  that  while  the  root-mean-square  (RMS) 
fluctuation  level  of  the  streamwiae  velocity  was  increased  by  the  presence  of  the 
polymer  solution,  the  RMS  of  the  wall-normal  velocity  and  the  Reynolds  shear  stress 
were  reduced  as  the  polymer  solution  mixed  with  the  flow.  Production  of  the 
streamwiae  Reynolds  normal  stress  was  decreased  but  the  production  of  the  Reynolds 
shear  stress  was  unchanged.  This  shows  that  the  pressure-strain  correlation  terms  in 
the  Reynolds  stress  transport  equations  may  be  directly  affected  by  the  polymers. 
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CHAPTER  1  -  INTRODUCTION 

The  phenomenon  of  drag  reduction  resulting  from  the  presence  of  long  chain, 
high  molecular  weight  polymers  in  turbulent  flow  of  liquids  has  received  much 
attention  in  the  past  several  decades.  Although  understanding  of  the  way  in  which 
polymers  reduce  drag  has  proved  elusive,  this  lack  of  understanding  has  not  impeded 
the  use  of  polymer  drag  reduction  for  engineering  purposes  in  pipe  flows  (see  the 
review  by  Hoyt,  1984).  One  area  receiving  little  attention  is  external  flows  where  the 
practical  implementation  of  polymer  drag  reduction  is  complicated  by  the  presence  of 
an  unbounded  fluid  medium.  In  these  flows,  the  use  of  polymers  requires  that  the 
polymer  solution  be  delivered  to  the  flow  and  the  performance  of  the  solution  as  a  drag 
reducer  will  depend,  to  large  extent,  on  the  mixing  of  the  polymer  with  the  ambient 
fluid. 

The  work  of  Wells  and  Spangler  (1967)  established  that  the  polymer  solution  had 
to  be  in  the  near-wall  region  to  cause  drag  reduction.  Investigations  by  Tiederman  et 
al.  (1985)  and  McComb  and  Rabie  (1982)  have  shown  that  drag  reduction  results  from 
the  presence  of  the  polymer  in  the  buffer  region  (10  <  y+  <  100).  Knowledge  of  the 
near-wall  structure  of  drag-reduced  turbulent  flows  is  confined  primarily  to 
homogeneous  flows  where  polymer  concentration  is  constant  and  the  turbulent  structure 
has  been  modified  throughout  the  flow.  In  flows  where  a  polymer  solution  is  injected  at 
the  wall,  one  can  examine  the  evolution  of  the  turbulent  structure  from  the  Newtonian 
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state  which  exists  upstream  of  the  injector,  to  a  drag-reduced  state  downstream.  In  this 
way,  the  effect  of  the  polymer  solution  on  the  near-wall  region  can  be  examined  in 
detail  and  die  initial  changes  to  the  near-wall  turbulent  structure  caused  by  the  polymer 
solution  can  be  detected. 

Although  the  evolution  of  the  mean  polymer  concentration  profile  downstream  of 
wall  injectors  has  been  investigated  by  various  authors,  the  mass  transport  processes 
that  govern  this  evolution  have  not  been  examined  directly.  The  Schmidt  number  for 
diffusion  of  polymer  solutions  in  water  is  high  (Sc  >  100)  hence,  mass  transfer  is 
dominated  by  turbulent  transport  mechanisms.  Since  turbulent  transport  is  convective 
in  nature  and  the  polymer  is  not  a  passive  scalar,  Le.  the  turbulence  is  modified  locally 
by  the  polymer,  mass  and  momentum  transport  in  injected  flows  are  coupled.  In  order 
to  model  this  type  of  flow,  direct  measurements  of  turbulent  mass  transport  and 
momentum  transport  are  necessary  unless  one  is  willing  to  proceed  with  adhoc 
assumptions. 

1.1  Objectives 

In  this  study,  a  turbulent  channel  flow  with  polymer  injection  at  the  wall  was 
examined.  The  specific  objectives  were  to  determine  in  detail  the  modifications  to  the 
structure  of  turbulence  caused  by  the  injection  of  the  polymer  solution  and  to  quantify 
the  high  Schmidt  number  mixing  process  which  occurs.  Realization  of  these  objectives 
required  two-component  velocity  measurements,  detailed  information  about  the 
polymer  concentration  field  and  direct  measurements  of  turbulent  mass  transport  The 
velocity  field  was  examined  using  two-component  laser  velocimetry.  To  determine  the 
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concentration  field,  a  laser-induced  fluorescence  technique  was  developed  to  measure 
instantaneous  polymer  concentration  profiles.  For  measurement  of  turbulent  mass 
transport,  the  concentration  measurements  were  combined  simultaneously  with  laser 
velocimetry  to  obtain  direct  measurements  of  turbulent  mass  transport 

In  what  follows,  the  behavior  of  the  polymer  solutions  is  discussed,  the  structure 
of  turbulence  in  drag-reduced  flows  is  outlined  and  the  relevant  work  in  mixing  and 
dispersion  for  flows  with  polymer  injection  is  examined. 

U  Drag-reducing  polymer  solutions 

Polymer  solutions  capable  of  causing  drag  reduction  in  turbulent  flows  are 
typically  non-Newtonian,  viscoelastic  fluids.  The  polymer  used  in  this  study  was  an 
aqueous  solution  of  SEP  ARAN  AP-273,  a  polyacrylamide  manufactured  by  Dow 
Chemical.  Argumedo  et  al.  (1978),  Tsai  and  Darby  (1978)  and  Cho  et  al.  (1983), 
among  others,  have  demonstrated  that  when  subjected  to  steady-shear  viscometric 
flows,  this  type  of  polymer  solution  exhibits  shear-thinning  behavior:  the  shear 
viscosity  decreases  with  increasing  shearing  strain-rate  (dU/dy).  Ryskin  (1987)  showed 
that,  when  placed  in  a  purely  elongational  strain  (dll/dx),  e.g.  flow  through  an 
axisymmetric  nozzle,  the  viscosity  increased  with  increasing  elongational  strain-rate. 
This  "extensional"  viscosity  can  be  up  to  two  orders  of  magnitude  larger  than  the  shear 
viscosity.  These  results  are  evidence  of  the  complex  constitutive  nature  of  polymer 
solutions.  Constitutive  relations  which  model  the  steady  shear  response  of  the  polymer 
solution  have  been  proposed  by  Argumedo  et  al.  and  Tsai  and  Darby.  Tsai  and  Darby 
also  modeled  the  response  of  the  fluid  to  certain  transient  inputs.  At  present  these 


rheological  modeling  efforts  are  at  best  tentative  and  have  not  been  implemented  in 
numerical  simulations  of  turbulent  flow. 
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In  flows  where  a  polymer  solution  is  injected,  the  rheology  of  the  resulting 
water-polymer  mixture  is  also  complicated  by  the  presence  of  concentration  variations. 
In  the  investigations  of  Frings  (1984)  and  Bewersdorff  (1984)  it  was  shown  that  when 
very- high-concentration  polymer  solutions  (Le.  several  thousand  parts  per  million)  are 
injected,  threads  of  high-concentration  fluid  can  form  which  persist  for  long  distances 
downstream.  This  nearly  "two-phase"  behavior  is  more  complex  ideologically  than 
flows  where  the  injected  concentration  is  below  1000  ppm  and  the  polymer  solution 
diffuses  more  readily.  To  date,  there  has  been  no  investigation  of  the  effect  of  non- 
homogeneous  concentration  on  the  ideological  properties  of  polymer  solutions. 


13  Turbulence  structure  and  drag  reduction 

For  Newtonain  flows,  the  turbulence  structure  is  defined  relative  to  the  Reynolds 
averaged  Navier-Stokes  and  continuity  equations.  For  a  stationary  turbulent  flow  these 
are: 
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The  terms  ujuj  are  unknowns  introduced  by  the  averaging  process.  These  unknowns 
comprise  the  Reynolds  stress  tensor  and  represent  transport  of  "i"  component 
momentum  in  the  "j"  direction  due  to  fluctuations  in  the  velocity  field.  'Turbulence 
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structure"  refers  to  the  spatial  variation  in  the  elements  of  the  Reynolds  stress  tensor 
and  other  statistical  quantities  appearing  in  Equation  1.1. 

In  flows  where  there  is  an  evolving  polymer  concentration  field,  the  right-hand 
side  of  Equation  1.1  could  contain  additional  terms  resulting  from  the  non-Newtonian 
nature  of  the  polymer  solution  and  the  unsteady  concentration  field.  Since  the  exact 
constitutive  equation  for  the  polymer  solution  has  not  been  determined,  the  form  of 
these  terms  is  unknown.  The  mathematical  form  of  the  Reynolds  stresses  (which  arise 
from  the  advective,  left-hand  side  terms  of  the  momentum  equation)  is  unaffected  by 
the  non-Newtonian  behavior  of  the  polymer  solutions.  Although  additional  terms  may 
be  introduced  when  polymer  is  present,  the  importance  of  die  Reynolds  stresses  is 

nnHiminisheri- 

Due  to  the  viscoelastic  nature  of  drag-reducing  polymers,  measurements  of 
velocity  using  hot-wire  anemometers,  hot-film  devices  or  pitot  tubes  can  be  grossly  in 
error.  This  leaves  the  laser  velocimeter  as  the  only  viable  method  for  measuring 
velocity  and  turbulence  statistics.  To  date  most  investigations  have  been  restricted  to 
one-component  velocity  measurements  although  recent  measurements  by  Luchik  and 
Tiederman  (1988),  and  Willmarth  et  al.  (1987)  have  quantified  drag-reduced  flows 
using  two-component  measurements.  These  measurements  have  been  confined  to  flows 
with  constant  polymer  concentration. 

The  turbulence  statistics  most  often  measured  are  profiles  of  the  mean  streamwise 
velocity  component.  These  have  been  reported  Berner  and  Scrivener  (1980),  McComb 
and  Rabie  (1982),  and  Bewersdorff  (1984)  for  circular  tubes,  Rudd  (1972)  and  Logan 


(1972)  for  square  tubes  and  Reischman  and  Tlederman  (1975),  Luchik  and  Tiederman 
(1988)  and  Willmarth  et  al.  (1987)  for  rectangular  channels.  The  effects  of  drag- 
reducing  additives  on  the  mean  velocity  profile  (normalized  with  inner  variables: 

U*  *  UAit,  y*  *  yut/v  where  u*  =  ‘sjx^/p  and  v  is  the  kinematic  viscosity  evaluated  at 
the  long  time  average  wall  shear  rate)  are  clearly  established.  The  linear  sublayer 
(y+  <  8)  remains  unaffected  and  the  logarithmic  overlap  region  is  displaced  toward 
higher  velocities  but  its  slope  is  unchanged  relative  to  Newtonian  flows. 

The  root  mean  square  (rms)  of  the  stream  wise  velocity  (uO  also  exhibits  specific 
trends  in  drag-reduced  flows.  In  Newtonian  flows  the  peak  in  uTu*  occurs  at  about 
y+  =  15;  for  drag  reducing  flows  this  peak  is  shifted  outward  in  both  physical  distance 
and  y+.  These  changes  can  be  directly  related  to  the  changes  which  occur  in  the  mean 
velocity  profile;  an  increase  in  the  mean  velocity  gradient  results  in  an  increase  in  u' 
and  the  opposite  also  holds. 

The  results  presented  by  Luchik  and  Tiederman  (1988)  and  Willmarth  et  al. 
(1987)  show  that  the  RMS  of  the  normal  velocity  decreases  both  physically  and 
nondimensionally  throughout  the  viscous  sublayer.  Similar  results  were  obtained  by 
Logan  (1972)  but  the  interpretation  of  those  results  was  complicated  by  the  presence  of 
secondary  flows  in  the  square  duct  used  in  that  study. 

Luchik  and  Tiederman  (1988)  also  measured  the  turbulent  shear  stress,  uv,  and 
found  that  the  location  of  the  peak  in  uv  moved  outward  and  the  peak  value  was 
decreased.  Willmarth  et  al.  (1987)  found  uv  greatly  decreased  across  the  entire  flow 
and  proposed  the  existence  of  a  non-Newtonian  retarding  force  to  account  for  the 
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discrepancy  between  the  total  stress  and  the  sum  of  the  turbulent  and  mean  viscous 
shear  stresses.  This  behavior  was  not  seen  in  the  experiments  of  Luchik  and 
Hederman,  where  the  polymer  was  mixed  uniformly  with  the  flow,  but  is  similar  to  that 
presented  by  Bewersdorff  (1984)  for  pipe  flows  where  a  very  high  concentration 
polymer  solution  was  injected  at  the  centerline.  In  Bewersdorff  s  flow,  the  solution  did 
not  completely  mix  with  the  ambient  fluid  resulting  in  the  presence  of  high- 
concentration  polymer  threads  in  the  flow.  This  raises  the  question  of  whether  polymer 
threads  were  present  in  the  flow  of  Willmarth  et  al. 

These  studies  indicate  that  the  presence  of  polymers  in  a  turbulent  flow  results  in 
significant  changes  in  the  structure  of  turbulence  which,  in  turn,  cause  reduced  viscous 
drag.  The  mechanism(s)  by  which  die  polymer  solutions  cause  these  structural  changes 
have  not  been  identified. 


1.4  Scalar  transport  and  the  dispersion  of  injected  polymer  solutions 


The  conservation  of  a  scalar  quantity  C  in  stationary  turbulent  flow  is  governed 
by  the  Reynolds  averaged  species  conservation  equation: 


a 

dxj 
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The  time-averaging  process  introduces  the  unknowns  ujc.  For  a  developing  scalar  field 
in  a  two-dimensional  (on  average)  boundary  layer  or  channel  flow  one  would  usually 
reduce  this  equation  to  its  "thin  shear  layer"  (or  boundary  layer)  form: 
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The  unknown  vc  represents  scalar  transport  in  the  direction  normal  to  the  wall  resulting 
from  velocity  and  concentration  fluctuations.  Solution  of  this  model  equation  requires 
the  specification  of  vc  as  a  function  of  position. 


Flows  with  injection  at  the  wall  occur  most  frequently  in  film-cooling 
applications  in  gas  turbines  and  other  harsh  environments.  The  review  by  Hartnett 
(1985)  shows  that  the  injection  flowrates  which  are  typical  for  film  cooling  are 
significantly  larger  than  those  which  are  necessary  for  polymer  drag  reduction.  This 
results  in  the  injected  fluid  extending  far  into  the  fully-turbulent  portion  of  the  flow. 
Since  this  radically  alters  the  near-wall  turbulent  structure,  the  physics  of  the  injection 
process  for  film  cooling  are  significantly  different  from  those  of  polymer  injection 
where  there  is  minimal  disturbance  to  the  flow. 


Other  studies  of  scalar  transport  in  Newtonian  flows  have  been  generally 
i  restricted  to  investigations  of  heat  transfer.  Hishida  et  al.  (1986)  examined  a  fully 

developed  pipe  flow  with  constant  heat  flux  at  the  wall.  That  study  showed  that  mean 
velocity  and  temperature  profiles  were  similar  and  turbulent  scalar  transport  (v0) 
profiles  were  similar  to  those  for  turbulent  momentum  transport  (uv).  Antonia  et  al. 
(1977)  examined  vfl  in  a  boundary  layer  at  several  locations  downstream  of  a  step 
change  in  wall  heat  flux.  Their  results  showed  that  the  maximum  value  of  v9  at  all 
locations  downstream  of  the  discontinuity  was  approximately  equal  to  the  heat  flux  at 
the  wall.  When  distance  from  the  wall  was  scaled  with  the  thermal  boundary  layer 

i 

thickness,  there  was  good  agreement  between  the  v6  profiles  for  all  streamwise 
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locations.  No  direct  comparisons  were  drawn  between  uv  and  v9;  however,  as  the 
thickness  of  the  thermal  boundary  layer  approached  that  of  the  momentum  boundary 
layer,  the  vd  profile  was  similar  to  typical  behavior  of  uv  in  boundary  layer  flows. 
These  results  were  shown  to  also  be  consistent  with  behavior  in  flows  with  a  step 
change  in  wall  temperature. 

The  invesigations  of  Hishida  et  al.  and  Antonia  et  al.  show  that  for  fully 
developed  flows  and  boundary  layer  flows  where  the  thermal  boundary  layer  thickness 
equals  the  momentum  boundary  layer  thickness,  turbulent  scalar  transport  profiles  are 
similar  to  turbulent  momentum  transport  profiles.  This  occurred  because,  in  both 
studies,  the  working  fluid  was  air  (where  the  Prandtl  number  is  near  unity)  and  the 
boundary  conditions  were  similar  for  the  momentum  and  energy  equations.  In  flows 
with  injection  at  the  wall,  this  similarity  does  not  exist:  The  boundary  condition  at  the 
wall  for  the  species  conservation  equation  specifies  zero  scalar  flux  and  for  mass 
transport  in  liquids,  the  Schmidt  number  is  at  least  several  hundred.  These  differences 
together  with  the  fact  that  the  flow  is  not  fully  developed  make  similarity  between 
turbulent  mass  transport  and  momentum  transport  unlikely. 

Investigations  into  the  mixing  of  injected  polymer  solutions  have  been  limited  to 
the  measurement  of  mean  concentration  profiles.  Latto  and  El  Reidy  (1976,1984) 
measured  mean  concentration  profiles  downstream  of  slot  injectors  in  a  flat  plate 
boundary  layer  and  showed  that  polymer  solutions  diffuse  more  slowly  than  water 
solutions.  An  investigation  by  Collins  and  Gorton  (1976)  showed  that  downstream  of 
the  region  where  large  concentration  gradients  exist,  diffusion  rates  are  comparable  to 
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those  resulting  from  plain  water  injection.  Fabula  and  Bums  (1970)  measured  the  mean 
polymer  concentration  profile  far  downstream  of  a  slot  injector  in  a  zero  pressure 
gradient  boundary  layer.  They  showed  that  in  this  region,  where  the  polymer  was 
present  across  the  entire  boundary  layer,  the  concentration  profile  was  self-similar. 

All  these  measurements  were  made  using  pitot  probes  and  isokinetic  sampling 
techniques.  Latto  et  al.  (1981)  have  shown  that  measured  mean  concentrations  are 
always  less  than  actual  concentrations  in  inhomogeneous  flows  when  this  measurement 
technique  is  used.  This  is  due  to  the  strong  dependence  of  solution  viscosity  upon 
concentration. 

None  of  these  polymer  studies  examined  turbulent  mass  transport  directly,  nor 
did  any  attempt  to  relate  the  evolution  of  the  concentration  field  to  the  changes  in  the 
structure  of  turbulence  which  result  in  drag  reduction. 

li  Closure 

In  the  next  chapter,  the  experimental  facilities  and  the  apparatus  and  procedures 
used  for  this  investigation  are  discussed.  This  includes  a  discussion  of  the  two- 
component  laser  velocimeter  system  used  for  making  high  spatial  resolution  velocity 
measurements  and  procedures  to  eliminate  bias  errors  which  can  occur  with  this 
technique.  A  laser-induced  fluorescence  technique  for  making  accurate,  high-speed, 
non-intrusive  measurements  of  scalar  concentration  profiles  is  described  in  detail.  And 
the  combination  of  laser  velocimetry  and  concentration  measurements  to  obtain  direct 
measurements  of  turbulent  mass  transport  is  discussed. 
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The  results  from  this  study  are  contained  in  the  next  four  chapters.  In  Chapter  3, 
measurements  of  turbulence  structure  for  the  water  flow  with  no  injection  are  presented. 
These  data  establish  the  "standard"  character  of  the  basic  water  flow  and  serve  as  a 
yardstick  by  which  modification  to  the  flow  can  be  measured.  Also,  a  highly  accurate 
method  for  determining  the  wall  shear  stress  in  fully  developed  flows  is  shown.  The 
concentration  field  in  injected  flows  is  discussed  in  Chapter  4.  Time-averaged  statistics 
are  presented  and  the  underlying  behavior  of  the  concentration  field  is  investigated  by 
examining  the  time-history  of  instantaneous  concentration  profiles.  A  comparison 
between  the  concentration  field  resulting  from  polymer  injection  and  that  resulting  from 
the  injection  of  dyed  water  is  made.  In  Chapter  5,  the  effect  of  the  evolving 
concentration  field  on  the  structure  of  turbulence  is  examined.  Changes  in  the  mean 
quantities  and  Reynolds  stresses  resulting  from  the  developing  concentration  field  are 
discussed,  and  the  causes  of  these  changes  are  inferred  from  examination  of  the 
Reynolds  stress  transport  equations.  In  these  results,  the  effect  of  the  injection  process 
on  turbulence  structure  is  isolated  from  the  modification  of  the  flow  by  the  polymers  by 
examining  a  flow  with  water  injection.  Chapter  6  presents  direct  measurements  of 
turbulent  mass  transport  for  both  water  injection  and  polymer  injection.  These  results 
are  compared  and  the  relationship  between  the  velocity  and  concentration  fluctuations 
is  examined  in  detail.  A  discussion  of  issues  relevant  to  modeling  these  flows  is  then 
presented. 

The  final  chapter  summarizes  the  major  results  and  presents  conclusions. 


CHAPTER  2  -  APPARATUS  AND  PROCEDURES 


2.1  Experimental  facilities 

The  water  flow  loop  used  in  these  experiments,  shown  in  Figure  2.1,  was  driven 
by  four  ninety  gallon  per  minute  centrifugal  pumps  operating  in  parallel  yielding  a 
maximum  Reynolds  number,  based  on  channel  height,  of  about  40  000.  At  each  end  of 
the  test  section  there  was  a  large  stilling  tank  to  isolate  the  test  section  from  any 
hydrodynamic  disturbances  in  the  flow  loop.  The  upstream  stilling  tank  contained  a 
perforated  plate  followed  by  a  screen  and  open-cell  sponge  section  and  a  smooth  two- 
dimensional  contraction  at  the  oudet.  The  inlet  of  the  channel  was  preceded  by  a  flow 
straightener  consisting  of  closely  packed  plastic  drinking  straws  which  insured  that  no 
large-scale  vorticity  existed  in  the  channel  entry  flow.  The  downstream  tank  contained 
a  perforated  plate  to  damp  out  disturbances  and  a  copper  coil  through  which  cooling 
water  was  passed  to  maintain  the  channel  water  at  a  constant  temperature. 

The  flow  loop  had  a  rectangular-cross-section  channel  (6.0  cm  high  by  57.5  cm 
wide)  as  the  test  section.  The  test  section  was  constructed  from  one- half  inch  acrylic 
and  polycarbonate  sheet  and  was  more  than  one  hundred  channel  heights  long. 
Polymer  solutions  were  injected  through  flush-mounted,  angled  slots  located  in  both  of 
the  57.5  cm  walls  of  the  channel.  The  injection  slots  had  a  width,  measured  in  the 
stream  wise  direction,  of  2.5  mm  and  were  inclined  25  degrees  to  the  flow  direction  as 
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shown  in  Figure  2.2.  The  injectors  were  located  about  eighty  channel  heights 
downstream  of  the  inlet 

The  coordinate  system  used  in  the  description  of  the  apparatus  and  presentation 
of  die  results  is  also  shown  in  Figure  22.  The  x-axis  is  parallel  to  the  flow  direction 
and  the  y-axis  is  normal  to  the  flow  direction  and  the  long  (57.5  cm)  dimension  of  the 
channel  cross  section.  The  origin  of  this  coordinate  system  is  located  at  the  center 
(mid-span)  of  the  injection  slot  located  in  the  lower  57.5  cm  wall  of  the  channel. 

22  Polymer  solution  preparation 

These  experiments  conducted  over  a  period  of  several  months  using  several 
different  batches  of  polymer  solution.  To  ensure  that  the  results  are  comparable,  the 
polymer  solutions  were  prepared  using  a  standard  procedure  and  the  batches  were 
checked  for  repeatability. 

The  additive  was  an  aqueous  solution  of  SEP  ARAN  ap-273,  a  polyacrylamide 
manufactured  by  Dow  Chemical,  at  a  concentration  of  700  ppm  based  on  weight.  The 
dry  powder  was  suspended  in  300  mO  of  isopropyl  alcohol  and  mixed  into  filtered, 
deaerated  tap  water  at  about  38  °C  using  a  magnetic  stirrer.  The  concentration  of  this 
initial  mixture  was  about  5000  ppm  and  it  was  allowed  to  hydrate  for  12  to  24  hours. 
This  solution  was  then  diluted  to  the  desired  concentration  using  filtered  tap  water  and 
allowed  to  hydrate  for  12  to  24  hours  before  using. 

The  repeatability  of  the  polymer  solutions  was  checked  in  two  ways.  The  steady 
shear  viscosity  was  checked  over  a  range  of  shear  rates  using  a  Brookfield  LVT-SCP3 
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Figure  2.2  Injection  slot  and  optical  arrangement 
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cone  and  plate  viscometer.  Also,  samples  of  each  batch  were  diluted  to  100  ppm  and 
drag  reduction  was  measured  in  a  1.405  cm  diameter  tube  for  a  range  of  flowrates. 
Only  polymer  solutions  which  exhibited  reasonable  consistency  in  these  two  measures 
were  used 

23  Velocity  measurements 

2.3.1  Two-component  laser  velocimeter  system 

The  transmitting  optics  consisted  of  a  Thermo-Systems  Incorporated  (TSI)  model 
9100-8  laser  velocimeter  system  incorporating  a  Lexel  model  85.5  500  mW  argon-ion 
laser.  This  is  a  two-component,  three-beam  system  where  velocities  are  measured 
using  one  blue  (488  nm)  beam,  one  green  (514.5  nm)  beam,  and  a  blue-green  beam 
containing  both  wavelengths.  Both  channels  were  frequency  shifted  40  MHz  using  a 
single  Biagg  cell  ahead  of  the  color  separator.  The  beam  spacing  was  25  mm  and  the 
beams  were  expanded  to  2.5  mm  and  focused  with  a  250  mm  transmitting  lens. 

Modifications  to  the  standard  TSI  9100-8  system  consisted  of  the  following:  The 
beam  displacer  in  the  output  end  of  the  model  9112-1  dichroic  color  separator  was 
rotated  90  degrees  to  place  the  green  beam  at  90  degrees  from  the  blue  beam.  After  the 
color  separator,  the  blue-green  beam  was  displaced  to  the  optical  axis  using  a  TSI  9174 
50  mm  beam  displacer.  These  changes  placed  the  blue-green  beam  on  the  optical  axis 
and  the  blue  beam  and  green  beam  at  nominally  ±  45  degrees  from  vertical.  In  this 
way,  velocity  components  at  about  ±  45  degrees  to  the  flow  direction  were  measured. 
With  this  system  the  included  angle  between  the  beams  can  deviate  from  90  degrees  by 
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as  much  as  4  degrees.  This  is  accounted  for  in  the  data  analysis  procedure  discussed 
below. 

For  this  system  the  measured  velocity  components  are  in  a  plane  which  is  parallel 
to  the  flow  direction  but  inclined  at  an  angle  of  about  3  degrees  relative  to  the  wall- 
normal  direction.  This  introduces  some  error  into  die  measurement  of  quantities  related 
to  the  normal  velocity.  For  the  second  moment  of  the  normal  velocity  v2, 

vm2  =  v2  cos2a  +  w2  sin2a  (2.1) 

where  vm2  is  the  second  moment  of  the  measured  velocity  component  and  a  is  the 
angle  of  inclination  from  the  wall-normal  direction.  For  the  turbulent  shear  stress  uv 

uvm  =  uv  cosa .  (2.2) 

Thus  the  error  introduced  in  uv  for  a  =  3  degrees  is  about  0.1  percent  which  is 
considerably  smaller  than  the  uncertainty  in  uvm.  The  effect  on  the  second  moment  of 
the  normal  velocity  can  be  evaluated  using  the  near-wall  values  of  v2  and  w2  from  Kim 

et  al.  (1987).  Table  2.1  shows  the  resulting  error  in  the  measured  value  of  v'  =  for 
a  =  3  degrees.  The  resulting  error  in  v'  is  only  significant  very  near  the  wall  (y+  <  5). 

The  optics  table  was  equipped  to  translate  the  measurement  volume  in  three 
orthogonal  directions.  The  beams  were  translated  in  the  vertical  (y)  direction  using  a 
pair  of  enhanced  aluminum  mirrors.  The  upper  mirror  and  the  transmitting  lens  were 
mounted  on  a  vertical  translation  stage  which  was  positioned  using  two  two-inch 
micrometers  (see  Luchik,  1985).  The  entire  vertical  traverse  assembly  was  mounted 
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Table  2.1  Estimated  error  in  wall-normal  root-mean-square  velocities  from  results  of 
Kim  etaL  (1987). 


y+ 

V  m 

v' 

%error 

1.344 

3.8158 

281.6 

1.934 

1.5876 

58.8 

2.630 

1.2829 

28.3 

3.433 

1.1517 

15.2 

4.341 

1.0432 

4.2 

5.354 

1.0287 

2.9 

6.472 

1.0209 

2.1 

7.639 

1.0152 

1.5 

9.018 

1.0113 

1.1 

10.440 

1.0116 

1.1 

on  a  translation  stage  which  moved  the  measurement  volume  parallel  to  the  optical 
axis  (z-dinection).  The  optics  table  was  located  on  translation  stages  which  allowed 
several  inches  of  travel  in  the  stream  wise  (x)  direction. 

The  location  of  the  measurement  volume  relative  to  the  wall  was  determined 
using  a  TSI  model  9140  receiving  optics  assembly  and  eyepiece  installed  in  the 
transmitting  optics  train.  By  moving  the  probe  volume  toward  the  wall,  the  y-location 
where  the  image  of  the  probe  volume  intersected  the  wall  could  be  determined  to  within 
about  one-half  the  probe  volume  diameter.  This  position  was  then  used  as  a  reference 
for  all  subsequent  vertical  location  measurements. 


The  beam  geometry  is  defined  by  the  transmitting  optics  and  determines  exactly 
which  components  of  velocity  are  measured  and  the  relationship  between  measured 
Doppler  frequency  and  velocity.  With  the  beams  in  air,  a  screen  was  placed  a  known 
distance  several  feet  from  the  beam  crossing  point  and  oriented  normal  to  the  optical 


i 
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axis.  The  points  where  the  three  beams  intersected  the  screen  were  marked  and  the 
lengths  of  the  sides  of  the  triangle  formed  by  these  points  was  measured.  From  these 
measurements  and  the  known  distance  from  the  beam  crossing  point,  the  half-angle  of 
each  beam  pair  was  calculated  and  used  to  determine  the  fringe  spacing  for  each 
channel  The  lengths  of  the  three  sides  of  this  triangle  are  used  to  calculate  the  included 
angle  between  the  measured  velocity  components  which  in  turn  is  used  to  calculate  true 
orthogonal  velocity  components  from  the  measured  velocities. 

The  receiving  optics  were  standard  TSI  two-component  items.  Scattered  light  was 
collected  using  a  250  mm  focal  length  lens.  A  2.27X  telescope  was  inserted  behind  the 
collection  lens  to  increase  the  spatial  resolution  of  the  receiving  optics.  A  50  pm 
aperture  was  used  in  a  TSI  9143  field-stop  assembly  to  ensure  spatial  coincidence  of  the 
two  measured  velocities.  The  light  passing  through  the  aperture  was  re-collimated  and 
passed  through  a  dichroic  color  separator.  The  separate  colors  were  bandpass  filtered 
and  focused  on  two  photomultiplier  tubes. 

The  spatial  resolution  of  the  measurement  system  is  defined  by  both  the 
transmitting  optics  and  the  receiving  optics  and  their  relative  orientation.  The  axis  of 
the  receiving  optics  was  located  in  the  y-z  plane  at  an  angle  of  9.5  degrees  from  the  axis 
of  the  transmitting  optics.  For  this  orientation,  the  extent  of  the  probe  volume  in  the  y- 
direction  is  defined  by  the  diameter  of  the  beam  intersection  (ly  =  63  pm).  The  x 
dimension  of  the  measurement  volume  is  determined  by  the  diameter  of  the  field  of 
view  of  the  receiving  optics  (lx  =  55  pm).  The  length  of  the  probe  volume  along  the  z- 
axis,  defined  by  the  field  of  view  of  the  receiving  optics  and  the  collection  angle,  was 
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lz  =  330  Jim.  Non-dimensionalizing  with  the  shear  velocity  and  kinematic  viscosity  of 
the  basic  water  flow  yields:  lx+  =  1.6,  ly+  =  1.8,  and  lz*  =  9.6. 

The  photomultiplier  output  for  each  channel  was  high-pass  filtered  at  25  MHz  to 
remove  die  pedestal  frequencies  and  was  electronically  down-mixed  to  yield  an 
effective  shift  frequency  of  0.5  MHz  using  a  TSI  model  9186  down-mixer.  The 
frequency  of  the  down-mixed  signals  was  determined  using  two  TSI  model  1980  counter 
processors.  Signals  were  high-limit  and  low-limit  filtered  at  1  MHz  and  30  kHz 
respectively  resulting  in  a  pass-band  equivalent  to  ±  2.5  m/s.  (The  center-line  velocity 
for  the  water  flow  was  0.65  m/s.)  The  Schmitt-trigger  threshold  in  these  counters  is 
50  mV  so  amplifier  gain  was  adjusted  to  maintain  noise  levels  below  10-15  mV. 
Maximum-signal  amplitudes  of  about  1  V  resulted  in  signal-to-noise  ratios  up  to  60:1. 
The  elapsed  time  for  8  zero  crossings  was  measured.  Validation  required  less  than 
three  percent  difference  in  the  period  determined  from  the  total  time  and  the  period 
determined  from  the  time  for  5  zero  crossings. 

Velocity  data  were  acquired  using  a  Masscomp  5520  microcomputer 
communicating  with  the  TSI  model  1998  interface  through  a  Masscomp  PI16F  fast 
parallel  interface.  Data  were  not  transferred  to  the  computer  unless  measurements  were 
made  on  both  channels  within  a  predetermined  coincidence  window.  The  length  of  the 
coincidence  window  was  set  equal  to  the  estimated  minimum  transit  time  through  the 
probe  volume  (about  100  jis). 
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Z3.2  Seeding,  sampling  and  velocity  bias 

The  major  drawback  to  measuring  fluid  velocities  using  laser  velocimetry  is  that 
the  fluid  velocity  is  not  measured  directly;  the  measured  velocity  is  that  of  particles 
suspended  in  die  fluid.  To  obtain  accurate  measurements  one  must  ensure  that  there  is 
no  more  than  one  particle  in  the  measurement  volume  at  any  given  time.  One  must  also 
ralce  care  to  ensure  that  the  estimated  velocity  statistics  are  not  affected  by  fringe  bias, 
filter  bias,  or  velocity  bias  (see  Edwards,  1987). 

For  this  study  the  flow  was  seeded  with  homogenized  whole  milk  which  contains 
fat  particles  about  0.3  }im  in  diameter.  To  insure  that  the  flow  was  uniformly  seeded, 
the  injected  fluids  were  seeded  at  the  same  concentration  as  the  water  in  the  flow  loop. 
In  the  case  of  the  polymer  solutions,  this  was  accomplished  by  adding  the  seed  to  the 
water  used  to  dilute  the  solution  to  its  final  concentration. 

The  maximum  particle  arrival  rate  achieveable  with  a  given  laser  velocimeter  in  a 
given  flow  (with  only  one  particle  in  the  measurement  volume  on  average)  can  be 
estimated  by  dividing  the  velocity  of  the  flow  by  the  diameter  of  the  probe  volume.  For 
the  flow  in  question  this  results  in  a  particle  arrival  rate  of  about  4  kHz  to  8  kHz, 
depending  on  position  in  the  flow.  For  the  results  presented  here,  the  water  in  the  flow 
loop  was  prepared  by  filtering  with  a  diatomaceous  earth  filter  to  remove  the  natural 
hydrosol  present  in  the  water.  The  water  was  then  seeded  with  whole  milk  at  a 
concentration  of  15  m0/1000  Q  which  resulted  in  maximum  validation  rates  of  about 


2  kHz.  This  insured  that  no  more  than  one  particle  was  in  the  probe  volume  at  any 
given  time. 
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Fringe  bias  was  overcome  by  frequency  shifting.  The  effective  fringe  velocity 
(after  down-mixing)  was  more  than  three  times  the  maximum  velocity  in  the  flow.  This 

I 

provided  almost  uniform  angular  response  for  the  velocity  measurements.  The 
possibility  of  filter  bias  was  eliminated  by  setting  the  filters  in  the  signal  processors  to 
yield  a  pass-band  equivalent  to  ±  2.5  m/s  (±  3  U0).  j 

Velocity  bias  occurs  because  there  is  a  higher  probability  of  measuring  a  high 
velocity  particle  than  a  low  velocity  particle  (see  McLaughlin  and  Tiederman,  1973).  j 

For  the  flow  in  question,  which  exhibits  thirty  to  forty  percent  turbulence  intensity  near 
the  wall,  the  error  in  the  estimate  of  the  mean  velocity  would  be  up  to  fifteen  percent  if 
this  effect  was  ignored.  For  the  two-component  measurements,  velocity  bias  was  I 

eliminated  using  the  fixed-waiting-time  sampling  method.  In  this  method  the  data 
acquisition  system  is  inhibited  for  a  fixed  amount  of  time  between  one  realization  and  ^ 

the  next  The  Kolmogorov  time  scale  is  representative  of  the  smallest  or  dissipation 
scales  in  the  flow.  (For  isotropic  turbulence,  this  time  scale  is  generally  smaller  than 
the  Eulerian  microscale.)  For  the  flow  in  question,  the  Kolmogorov  time  scale  was  I 

estimated  by  assuming  that  dissipation  of  turbulence  kinetic  energy  was  equal  to 
production  at  the  production  peak.  The  validation  rate  was  maintained  at  about  five 

I 

times  the  Kolmogorov  frequency  and  the  inhibit  time  was  ten  times  the  inverse  of  the 

validation  rate.  Based  on  the  results  of  Gould  et  al.  (1988),  this  approach  should  be 

adequate  to  eliminate  velocity  bias  effects  from  the  mean  quantities  and  it  is  assumed  I 

that  the  effects  on  higher  statistical  moments  are  also  eliminated. 


23 


For  the  water  flow  results  presented  in  Chapter  3,  all  statistics  were  based  on 
ensembles  of  25  000  data  realizations.  Statistics  for  the  flows  with  water  and  polymer 
injection  (Chapter  5)  were  calculated  using  ensembles  of  10  000  realizations. 

1A  Concentration  measurements 

Time-resolved  concentration  measurements  were  obtained  in  a  manner  similar  to 
that  used  by  Koochesfahani  and  Dimotaltis  (1986).  The  injected  fluid  was  marked  with 
a  fluorescent  dye  and  the  spatial  distribution  of  the  intensity  of  fluoresced  light  emitted 
from  a  laser  beam  directed  normal  to  the  channel  wall  was  measured.  The  dye 
concentration  at  a  point  was  determined  from  the  fluorescent  light  intensity  and  the 
injectant  concentration  was  inferred  from  the  measured  dye  concentration.  For 
turbulent  flows,  where  dispersion  occurs  due  to  convective  mixing,  the  "turbulent"  mass 
diffusivity  is  typically  two  or  more  orders  of  magnitude  larger  than  the  molecular  mass 
diffusivity.  This  indicates  that  the  time  scale  for  molecular  diffusion  is  more  than  one 
hundred  times  the  timescale  for  "turbulent  diffusion",  hence  the  effect  of  molecular 
diffusion  will  be  small  and  the  dye  concentration  will  yield  a  good  estimate  of  the 
polymer  concentration. 

2.4.1  Analysis 

For  light  propagating  in  an  absorbing  medium,  the  change  in  intensity,  die,  over 
the  differential  distance  dy  is 


die  =  -kle  dy 


(2.3) 
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where  k  is  the  extinction  coefficient  for  the  medium.  In  this  case,  Ie  represents  the 
intensity  of  the  of  the  laser  beam  used  to  excite  the  dye.  Integration  of  Equation  2.3 
from  the  point  where  the  beam  enters  the  medium  (y  *  0)  to  the  point  in  question  (y) 
yields 

J 

-  fk(r)  dr 

i«(y)  =  io  e  4  (2.4) 

where  Iq  =  I*(0).  The  intensity  of  the  fluoresced  light,  I,  emitted  from  a  differential 
length,  ds,  along  the  excitation  beam  is  some  fraction  of  the  light  absorbed  over  that 
distance, 

I  =  adle=akleds  (2.5) 

where  a  is  a  constant  Hence,  for  a  small  segment  of  the  excitation  laser  beam  of 
length  d  at  location  y, 

I(y)  =  ak(y)Ie(y)d .  (2.6) 

For  fluorescent  dyes  in  the  range  of  concentrations  of  interest  here,  the  extinction 
coefficient  is  proportional  to  the  dye  concentration 

k  =  (3  C  .  (2.7) 

Combining  the  above  equations  yields 

7 

-fPC(r)dr 

I(y)  =  a  J3C(y)  4  e  “  d.  (2.8) 

This  equation  relates  the  fluorescence  intensity  from  a  discreet  segment  of  the  laser 
beam  to  the  initial  laser  beam  intensity  and  the  dye  concentration  along  the  beam  path. 
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Solving  for  the  local  concentration  C(y)  yields 


C(y)  = 


I(y) 

Iotxpd 


(2.9) 


If  the  dye  concentration  is  considered  constant  over  a  distance  s  measured  from 
the  midpoint  between  segments  i-1  and  i  to  the  midpoint  between  segments  i  and  i+1, 
then  the  integration  in  the  above  equation  can  be  written  as  a  summation.  The  average 
dye  concentration  over  the  segment  located  at  yn,  (the  nA  segment  of  the  laser  beam),  is 
given  by 


In 

Cn  =  C(yn)  =  -r— e  w 


Ioapd 


where  In  =  I(yn).  This  result  can  be  rewritten  as 


Cn  = 


IoCxPd 


ne^s 

F=1 


(2.10) 


(2.11) 


This  last  equation  expresses  the  dye  concentration  at  the  n*  location  along  the  beam 
path  in  terms  of  the  properties  of  the  dye,  the  dye  concentration  at  the  previous  n-1 
locations,  and  the  intensity  of  the  fluorescent  emission  from  a  small  segment  of  the 
excitation  beam.  Also  required  are  the  values  of  d  and  s  which  are  determined  by  the 
detector  used  for  the  fluorescence  intensity  measurements. 


For  these  experiments  the  fluorescence  intensity  distribution  along  the  path  of  the 
excitation  laser  beam  was  measured  using  a  Reticon  camera  incorporating  a  line  array 
of  photodiodes.  In  developing  working  equations,  the  characteristics  of  the  photodiode 
array  and  the  imaging  optics  must  be  considered.  The  voltage  for  a  given  element  of 


I 

the  array  increases  linearly  with  incident  light  intensity  from  a  dark-level  voltage.  Both 
the  dark-level  voltage  and  the  sensitivity  vary  slightly  from  element  to  element.  In 
addition,  since  measurements  were  made  near  a  wall,  the  fraction  of  the  total 
fluorescent  emission  collected  by  the  camera  lens  varied  with  position  in  the  field  of 
view. 


The  voltage  for  an  element  of  the  photodiode  array  is  related  to  the  fluorescent 
intensity  at  the  corresponding  point  in  the  field  of  view  by 


*  Bnhn  In  + 


(2.12) 


where  an  is  the  fraction  of  the  emitted  light  which  is  incident  on  the  element,  b„  is  the 
sensitivity  of  the  element  and  Vjn  is  the  dark  voltage  for  the  element.  Combining  this 
result  with  Equation  2.1 1  yields 


Q  = 


1 

Iotxpdanbn 


n— 1 

ne 

i=l 


pCjS 


(Vn  -  V^) 


(2.13) 


which  relates  concentration  to  measured  voltage.  Examination  of  this  equation  reveals 
that  the  term  in  brackets  depends  entirely  on  known  constants  and  element  number. 
Hence,  this  term  can  be  replaced  by  a  single  variable  An,  resulting  in 


n-1 

ne1 

i=l 


pciS 


(Vn-Vdn). 


(2.14) 


The  single  element-dependent  constant  is  then  determined  by  placing  a  known  dye 
concentration  in  the  flow  field  and  measuring  the  array  element  voltages. 
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2.4.2  Optical  arrangement 

For  these  experiments  the  fluorescence  intensity  distribution  along  the  path  of  the  j 

excitation  laser  beam  was  measured  using  a  Reticon  line-scan  camera  incorporating  a 
256  element  linear  photodiode  array.  The  elements  of  the  photodiode  array  are  25  pm 
by  425  pm  and  are  located  on  50  pm  centers.  For  all  stream  wise  measurement  stations  ^ 

in  the  flow  with  water  injection  and  the  first  measurement  station  in  the  flow  with 
polymer  injection,  the  magnification  was  set  to  yield  effective  measurement  volumes  . 

which  were  spaced  60.5  pm  apart  and  measured  30.3  pm  in  the  direction  normal  to  the 
wall.  The  magnification  at  the  other  stream  wise  locations  for  the  polymer  flow  resulted 
in  35.0  pm  measurement  volumes  spaced  70.0  pm  apart.  I 

The  blue  (488  nm)  line  from  a  Lexel  model  85.5  argon-ion  laser  was  used  to 
excite  the  dye.  The  diameter  of  the  excitation  beam  determined  the  spatial  resolution  of  I 

the  measurement  in  the  streamwise  and  spanwise  directions.  Focusing  the  excitation 
beam  with  a  200  mm  focal  length  lens  yielded  a  beam  diameter  of  115  ±  5  pm  across 
the  field  of  view.  1 

! 

f 

2.4.3  Data  acquisition 

I 

■  All  timing  and  data  acquisition  tasks  were  accomplished  using  a  Masscomp  5520 

: 

i 

j  micro-computer.  This  computer  has  8  Mbytes  of  random  access  memory  available  for 

I 

|  data  acquisition  and  sufficient  hardware  floating  point  calculation  capability  and  I 

I  graphics  capability  to  analyze  data  on-line  during  the  course  of  an  experiment.  Data 

| 

i  were  acquired  using  a  12  bit  A/D  converter  capable  of  one  million  conversions  per 
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signals  were  generated  using  the  clock  module  of  the  computer. 

The  Reticon  line-scan  camera  requires  a  start  signal  to  begin  each  scan  of  the 
array  and  a  clock  signal  to  drive  the  multiplexer  which  samples  die  individual  elements. 
The  output  from  each  element  is  in  die  form  of  a  dark-level  voltage  plus  a  voltage 
proportional  to  the  integrated  intensity  of  the  light  incident  on  the  element  during  the 
time  between  successive  scans  of  the  array.  Hence,  die  exposure  time  for  each  intensity 
profile  is  equal  to  the  time  between  start  signals.  The  maximum  output  voltage  is  two 
volts  and  the  dark  level  voltage  is  typically  about  10  mV. 

For  all  measurements  presented  here  the  camera  clock  signal  frequency  was  set  to 
857  kHz.  During  water  injection,  the  array  scan  rate  (the  frequency  at  which  start 
signals  woe  delivered)  was  3  kHz  yielding  an  exposure  time  of  333  microseconds. 
Every  sixth  scan  of  the  array  was  digitized  resulting  in  a  sampling  rate  of  500  Hz.  For 
the  flow  with  polymer  injection,  the  array  scan  rate  was  2  kHz  yielding  an  exposure 
time  of  500  microseconds.  Two  sampling  rates  were  used  in  the  polymer  injection 
flow.  To  calculate  long-time  averaged  statistics,  every  tenth  scan  of  the  array  was 
digitized  yielding  a  sampling  rate  of  200  Hz.  For  time-resolved  measurements,  every 
second  scan  of  the  array  was  digitized  for  a  sampling  rate  of  1000  Hz. 

2.4.4  Concentration  measurement  procedure 

Prior  to  the  concentration  measurements  the  injected  fluid  was  dyed  with 
fluorescein  disodium  salt  at  a  molar  concentration  of  about  1( T6.  The  water  in  the  flow 
loop  was  given  a  dye  concentration  of  about  3xl0"7  M  for  purposes  of  calibration. 
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The  fluorescence  of  the  fluorescein  dye  is  very  sensitive  to  pH  variations  (see  Walker, 

1986).  For  the  water  injection  measurements,  this  posed  no  problem  since  the  injected 
fluid  was  the  same  as  the  fluid  in  the  channel.  However,  the  polymer  solutions  had  a 
pH  of  8.1  whereas  the  channel  water  had  a  pH  of  7.4.  To  obtain  accurate  polymer 
concentration  measurements  the  pH  of  the  channel  water  was  adjusted  to  match  the 
polymer  solution  by  adding  sodium  hydroxide  to  the  water. 

The  fluorescent  emission  from  the  dye  in  the  polymer  solution  was  observed  to 
decrease  steadily  with  time.  Since  this  behavior  was  not  seen  when  water  was  the 
solvent,  it  is  believed  that  the  polymers  promoted  the  formation  of  non-fluorescing 
aggregates  of  dye  molecules.  Gentle  stirring  of  the  dyed  polymer  solution  slowed  this 
process  enough  that  there  was  negligible  change  in  fluorescence  during  the  course  of  an 
experiment 

The  camera  was  positioned  and  aligned  visually  with  the  excitation  beam  using 
the  through-the-lens  viewfinder.  Final  alignment  was  accomplished  by  observing  the 
camera  output  on  an  oscilloscope  and  adjusting  the  vertical  camera  position  to  place  the  I 

maximum  measured  light  level,  which  occurs  where  the  beam  first  enters  the  flowfield, 
at  the  first  photodiode. 

All  dye  concentrations  were  determined  using  an  attenuation  cell.  This  cell 
allowed  the  passage  of  a  laser  beam  through  a  cavity  in  which  a  fluid,  either  the 
injectant  or  water  from  the  flow  loop,  was  flowing.  By  measuring  the  attenuation  of  a 
laser  beam  passing  through  the  cell,  the  dye  concentration  of  the  fluid  could  be 
accurately  determined.  The  attenuation  measurement  was  performed  using  a  Lexel 
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model  504  laser  power  meter.  In  this  way,  the  dye  concentration  of  the  injectant  was 
measured  prior  to  an  experiment  and  the  dye  concentration  of  the  channel  water  was 
monitored  during  the  experiment 

To  avoid  buildup  of  dye  in  the  flow  loop,  dyed  fluid  was  injected  through  only  the 
center  60%  of  the  lower  injection  slot  and  undyed  fluid  was  injected  through  the  outer 
20%  on  each  end.  Undyed  fluid  was  injected  through  the  slot  in  the  top  channel  wall. 
This  combination  resulted  in  the  average  dye  concentration  of  the  injected  fluid  being 
nearly  equal  to  the  background  concentration  in  the  channel  and  thus  kept  the 
background  concentration  relatively  constant 

Since  there  was  a  small  amount  of  fluctuation  in  the  laser  light  intensity  level 
(about  one  percent  at  a  frequency  of  about  50  Hz),  the  beam  was  sampled  before 
entering  the  test  section  and  a  photodiode  was  used  to  monitor  the  laser  intensity  during 
the  fluorescence  measurements.  The  signal  from  this  photodiode  was  sampled  and 
stored  with  each  scan  of  the  photodiode  array  in  the  camera.  This  information  was  used 
to  correct  all  fluorescence  measurements  for  fluctuations  in  the  excitation  source 
intensity. 

Three  separate  measurements  with  the  photodiode  array  were  required  for  each 
data  set  The  first  was  a  dark-level  measurement  in  which  1000  scans  of  the  array  were 
digitized  with  the  lens  aperture  covered  and  then  the  average  dark-level  voltage  for 
each  element  was  calculated  and  stored.  Simultaneously,  the  dark  level  for  the  beam- 
sampler  photodiode  was  measured  and  the  information  was  stored.  Next,  the 
fluorescence  intensity  resulting  from  the  presence  of  dye  in  the  channel  water  was 


k 


I 


1 


measured  with  the  camera  and  the  dye  concentration  of  the  water  in  the  flow  loop  was 
determined  using  the  attenuation  cell.  These  measurement  allowed  calculation  of  the 
constants  relating  the  voltage  from  each  element  of  the  photodiode  array  to  the  dye 
concentration  in  the  field  of  view.  Finally,  the  fluorescence  intensity  was  measured 
during  injection  with  die  photodiode  array.  These  steps  were  repeated  at  each 
streamwise  location  and  for  the  different  fluids  injected. 

All  concentration  statistics  were  calculated  using  ensembles  of  10  000 
realizations. 

2.4.5  Verification 

To  verify  the  accuracy  of  the  concentration  measurement  technique,  dyed  water 
was  injected  at  a  flow  rate  large  enough  to  establish  a  concentration  boundary  layer 
about  5  mm  thick  along  one  wall  of  the  channel.  The  intensity  of  the  attenuated  laser 
beam  leaving  the  test  section  was  measured  along  with  the  fluorescence  intensity 
profiles  along  the  beam.  The  average  concentration  along  the  beam  path  could  then  be 
deduced  from  the  measured  attenuation  of  the  excitation  laser  beam  and  from 
integration  of  the  measured  concentration  profile.  The  correlation  coefficient  for  these 
two  independent  measurements  of  space- averaged  concentration  was  0.97  verifying  the 
accuracy  of  the  concentration  profile  measurements.  A  portion  of  the  time  history  of 
these  two  measurements  is  shown  in  Figure  2.3. 
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24  Simultaneous  velocity  and  concentration  measurements 

Simultaneous  measurements  of  normal  velocity  and  concentration  at  a  point  were 
accomplished  by  combining  die  laser  induced  fluorescence  (L1F)  concentration 
measurement  technique  described  above  with  a  single-component  laser  velocimeter. 
When  combining  these  two  techniques,  two  considerations  are  paramount:  the  beams 
for  the  laser  velocimeter  must  not  affect  the  concentration  measurements;  and  the 
scattered  light  far  the  velocity  measurements  must  not  be  overwhelmed  by  the 
fluorescent  emission  from  the  dye.  In  practice,  these  goals  were  accomplished  by  using 
a  helium-neon  laser  for  the  velocity  measurements.  The  dye  used  for  the  LIF 
measurements  does  not  absorb  light  in  the  red  portion  of  the  spectrum,  hence  the 
helium-neon  laser  did  not  excite  the  dye  and  contaminate  the  concentration 
measurements.  The  fluorescent  emission  from  the  dye  is  primarily  in  the  green-yellow 
portion  of  the  spectrum,  away  from  the  helium-neon  wavelength  and  the  two  signals 
were  separated  effectively  using  an  optical  filter. 

2.5. 1  One-component  laser  velocimeter  system 

The  transmitting  optics  for  the  one-component  laser  velocimeter  consisted  of 
optics  from  the  previously  described  TSI  model  9100-8  system  used  in  conjunction  with 
a  Spectra-Physics  Inc.  model  106-1  10  mW  helium-neon  laser.  The  optical 
arrangement  used  is  shown  schematically  in  Figure  2.3.  A  TSI  model  9115-1  beam 
splitter,  set  for  reference-beam  operation,  was  used  to  create  two  equal-intensity, 
orthogonally  polarized  beams.  The  plane  of  the  beams  was  oriented  to  measure  the 
vertical  (y-component)  velocity.  The  polarization  of  one  beam  was  rotated  using  a 
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quarter-wave  reorder  to  match  that  of  the  other.  One  beam  was  frequency  shifted 
40  MHz  using  a  TSI  model  9182-12  Bragg  cell.  The  lower  beam  of  the  pair  was 
displaced  to  the  optical  axis  using  a  TSI  9174  beam  displacer,  to  enable  measurements 
in  close  proximity  to  the  lower  wall  of  the  channel.  The  beam  spacing  was  25  mm  and 
the  beams  were  expanded  to  1.54  mm  and  focused  with  a  250  mm  transmitting  lens.  As 
with  the  two-component  system  the  measured  velocity  component  is  inclined  at  about 
3  degrees  to  the  wall-normal  direction  and  die  errors  introduced  are  similarly  small. 

The  receiving  optics  were  the  same  as  the  two-component  system  except  they 
incorporated  only  one  photomultiplier  tube  and  a  narrow  bandpass  optical  filter 
centered  at  632.8  nm.  The  spatial  resolution  of  the  measurement  system  was 
unchanged  in  the  x-  and  z-directions,  but  die  smaller  beam  diameter  and  longer 
wavelength  of  the  helium-neon  system  resulted  in  a  vertical  probe  volume  dimension 
(ly)of  130pm(ly+  =  3.7). 

The  photomultiplier  output  was  high-pass  filtered  at  25  MHz  to  remove  the 
pedestal  frequencies  and  was  electronically  down-mixed  to  yield  an  effective  shift 
frequency  of  50  kHz  using  a  TSI  model  9186  down-mixer.  The  frequency  of  the  down- 
mixed  signals  was  determined  using  a  TSI  model  1980  counter-type  signal  processor. 
Signals  were  high-limit  and  low-limit  filtered  at  100  kHz  and  3  kHz  respectively 
resulting  in  a  pass-band  equivalent  to  ±0.3  m/s.  (The  maximum  RMS  normal  velocity 
was  about  0.03  m/s.)  Signal  quality  and  validation  requirements  were  the  same  as 
outlined  for  the  two-component  system.  An  analog  velocity  trace  was  generated  using 
a  TSI  model  1988  digital  to  analog  converter. 
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Seeding  was  the  same  for  this  system  as  for  the  two-component  system  and  the 
potential  sources  of  bias  were  either  eliminated  or  minimized.  Fringe  bias  was 
eliminated  by  frequency  shifting.  The  resulting  fringe  velocity  was  more  than  ten  times 
the  maximum  normal-velocity  component  and  it  was  large  enough  to  obtain 
measurements  on  the  particles  traveling  in  the  stream  wise  direction  at  the  maximum 
stream  wise  velocity.  Velocity  bias  was  again  eliminated  by  controlled  sampling  of  the 
processor  output  Instead  of  fixed-waiting-time  sampling,  as  used  for  the  two- 
component  measurements,  the  analog  velocity  signal  was  sampled  at  equal  time 
intervals.  As  before,  the  validation  rate  was  about  five  times  the  Kolmogorov 
frequency  and  the  sampling  rate  was  set  at  one-tenth  the  validation  rate.  According  to 
Edwards  (1987)  this  is  adequate  to  eliminate  Vv  incity  bias. 

2.5.2  Experimental  procedure 

In  performing  this  measurement,  two  conditions  had  to  be  satisfied.  The 
measurements  of  concentration  and  velocity  must  be  made  at  the  same  point  in  space; 
they  must  also  occur  at  essentially  the  same  point  in  time. 

The  procedure  used  to  ensure  spatial  coincidence  of  the  velocity  and 
concentration  measurements  consisted  of  the  following  steps.  Using  a  TSI  model  9140 
receiving  optics  assembly  and  eyepiece  installed  in  the  transmitting  optics  train,  an 
image  of  the  probe  volume  was  focused  on  the  eyepiece  reticule  and  aligned  with  the 
reticule  cross-hairs.  The  excitation  beam  for  the  concentration  measurements  was  then 
adjusted  so  that  it  remained  sharply  focused  on  the  reticule  and  centered  on  the  cross¬ 
hairs  throughout  the  vertical  travel  of  the  transmitting  optics.  The  line-scan  camera  was 
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then  aligned  and  focused  on  the  excitation  beam.  The  wall  location  was  established 
with  the  laser  velocimeter  using  the  procedure  outlined  for  the  two-component  system. 
The  exact  vertical  position  of  the  concentration  measurement  volumes  was  determined 
by  increasing  the  integration  time  of  the  video  camera  until  the  laser  velocimeter  beam 
crossing  registered  a  measurable  voltage  on  the  camera  output  Determining  the 
camera  element  numbers  which  corresponded  to  two  known  positions  of  the  beam 
crossing  accurately  established  the  position  of  the  concentration  measurement  volumes 
relative  to  the  channel  wall. 

The  velocity  signal  was  an  additional  analog  channel  for  the  A/D  converter  and 
the  velocity  was  sampled  once  for  each  intensity  profile  digitized.  A  block  diagram  of 
the  electronics  is  shown  in  Figure  2.4.  The  integration  time  for  the  concentration 
measurements  was  500  }is  which  is  one-fifth  the  estimated  smallest  timescale  in  the 
flow.  This  also  corresponds  to  the  average  time  between  validations  for  the  velocity 
measurements.  This  indicates  that  on  average  there  was  one  velocity  realization  per 
concentration  profile  measurement  thus  ensuring  good  temporal  coincidence.  Data 
were  acquired  at  a  rate  of  200  Hz  by  digitizing  every  tenth  concentration  profile 
measurement  and  ensembles  of  10  000  data  points  were  used  to  calculate  statistics. 

2.6  Data  analysis 

The  calculation  of  statistical  quantities  was  the  primary  data  analysis  task  for 
most  of  the  quantities  measured  in  this  study.  For  most  of  these  quantities,  this  process 
is  straightforward.  However,  the  two-component  laser  velocimeter  data  required 
somewhat  more  processing  and  warrants  a  more  detailed  discussion.  The  analysis  of 


Figure  2.5  Block  diagram  of  electronics  for  simultaneous  velocity  and  concentration 
measurement. 
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these  velocity  data  can  be  described  as  several  discrete  steps,  although  in  practice  some 
of  these  were  done  concurrently  to  reduce  computation  time. 

The  individual  velocities,  U'i  and  U'2,  measured  on  each  channel  were 
calculated.  Each  U'i.U^  pair  was  resolved  into  precisely  orthogonal  velocity 
components  Ui  and  U2  as  follows: 


Ui  =  v\ 

(2.15) 

u2  =  A  u'2-  *  uv 

sin0j  tan0j 

(2.16) 

In  these  equations,  6i  is  the  angle  between  the  measured  velocity  components 
determined  from  the  beam  geometry  measurements. 

All  statistics  were  calculated  using  the  orthogonal  U[,  U2  components.  The 
orientation  of  these  velocity  components  relative  to  the  flow  direction  was  determined 
from  the  velocity  measurements  themselves.  For  any  given  0  measured  relative  to  the 
flow  direction 


U  =  U2COS0  +  Ui  sin0 
V  =  U2Sin0  -  UicosO. 
Time-averaging  Equation  2.18  and  setting  V  =  0  yields 

0  =  Tan'1 


U2 


(2.17) 

(2.18) 


(2.19) 


The  value  of  0  was  determined  from  measurements  in  the  outer  flow  (where  maximum 
signal  quality  was  attained)  and  this  value  was  used  at  all  y-locations.  The  use  of  one 
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value  at  all  y-locations  resulted  in  average  normal-velocities  which  were  equal  to  zero 
within  the  uncertainty  of  the  normal-velocity  measurement. 

Equations  2.17  and  2.18  can  be  used  to  calculate  values  of  U  and  V  in  the  rotated 
coordinate  system.  Calculation  of  the  second  moments  in  the  rotated  coordinate  system 
is  more  complicated.  The  appropriate  relationships  are: 


u2  =  ui2sin20  +  U22cos20  +  2uiu2cos8sin0 

(2.20) 

v2  =  ui2cos20  +  U22sin20  -  2uju2cos0sin0 

(2.21) 

:  ^U22  -  Ui2  j  sin0cos0  +  ulul  ^sin20  -  cos20j 

(2.22) 

The  derivation  of  these  relationships  are  given  in  Appendix  A. 

The  calculation  of  some  quantities  required  either  integration  of  the  statistical 
quantities  across  the  channel  half-height  or  differentiation  of  quantities  at  a  point  in  the 
flow.  Integrations  were  done  using  the  trapezoidal  rule  and  assuming  zero  velocity  at 
the  wall.  Since  the  spaces  between  the  locations  of  the  measured  velocities  varied  with 
distance  from  the  wall,  derivatives  were  defined  as  the  average  of  the  one-sided  finite 
difference  approximations  on  either  side  of  the  point  in  question.  At  the  first  and  last 
data  point  in  a  profile,  the  derivatives  were  set  equal  to  the  one-sided  finite  difference 
approximations. 

The  determination  of  wall  shear  stress  can  be  done  in  several  different  ways.  For 
fully  developed  duct  flows,  the  average  wall  shear  stress  is  proportional  to  the  pressure 
gradient.  If  the  wall  shear  stress  is  constant  around  the  perimeter  of  the  duct  (as  it  is  in 
a  circular  pipe)  then  this  method  is  uncomplicated.  In  rectangular  channels,  the  wall 
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shear  stress  is  not  constant  around  the  perimeter  and  locally  the  wall  shear  stress  can 
deviate  from  the  average.  The  local  wall  shear  stress  can  be  determined  from  mean 
velocity  measurements  in  two  ways:  Estimating  the  mean  velocity  gradient  at  the  wall 
or  estimating  the  slope  of  the  logarithmic  portion  of  the  mean  velocity  profile.  Both  of 
these  methods  use  only  a  limited  amount  of  data  in  a  small  portion  of  the  flow. 

In  these  experiments,  the  shear  velocity  u*  =  ^x^/p  for  the  water  flow  with  no 
injection  was  determined  using  the  fact  that  the  total  shear  stress  (the  sum  of  the  mean 
viscous  and  turbulent  shear  stress)  must  be  a  linear  function  of  distance  from  the  wall 
for  a  fully  developed  flow: 


ttoui  =  U  4^-  -  puv  =  Tw»u(l-y/a) 
dy 


(2.23) 


where  a  is  the  channel  half-height.  Rearranging  yields 


which  expresses  a*  as  a  function  of  distance  from  the  wall.  The  shear  velocity  (which 
is  a  constant)  at  each  y-location  is  determined  using  the  measured  values  of  uv  and  the 
estimated  values  of  the  mean  velocity  gradient  for  each  measurement  point  in  the  flow. 
These  values  are  then  averaged  to  determine  the  shear  velocity  for  the  flow. 
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In  the  following  chapters,  experimental  data  will  be  presented.  For  mean  and 
and  root-mean-square  statistics,  the  uncertainty  in  the  estimated  quantity  is  comparable 
to  the  size  of  the  plotting  symbol  used.  For  statistics  such  as  uv  and  vc,  the  uncertainty 
estimates  are  indicated  by  error  bands  on  a  few  representative  points.  These  estimates 
are  at  the  95  percent  confidence  level  and  their  determination  is  discussed  in  Appendix 
B. 


CHAPTER  3  -  WATER  FLOW  RESULTS 


Measurements  in  the  water  flow  with  no  injection  were  made  100  mm 
downstream  of  the  injector.  These  results  are  presented  to  verify  the  "standard"  nature 
of  the  water  flow  and  as  a  reference  for  measurements  made  during  injection.  This  type 
of  flow  has  been  well  studied  by  many  investigators.  Hussain  and  Reynolds  (197S) 
performed  one -component  hot-wire  velocity  measurements  in  an  air  flow  in  a  32:1 
aspect  ratio  channel  which  minimized  end- wall  effects.  Kreplin  and  Eckelmann  (1979) 
examined  the  near-wall  structure  of  a  oil  channel  flow  where,  due  to  the  high  viscosity, 
excellent  spatial  resolution  was  obtained  with  hot  film  probes.  The  thesis  of  Wei  (1987) 
shows  high  spatial  resolution,  two-component  laser  velocimeter  measurements  over  a 
range  of  Reynolds  numbers.  Recently,  Kim  et  al.  (1987)  presented  results  from  a  full 
numerical  simulation  of  turbulent  channel  flow. 

3.1  Experimental  conditions 

Integral  parameters  were  calculated  and  are  presented  in  Table  3.1  along  with  the 
values  of  the  major  flow  variables.  Integration  of  the  mean  velocity  profile  yielded  a 
mass-averaged  velocity  of  0.574  m/s.  The  Reynolds  number  based  on  mass-averaged 
velocity  and  channel  height  was  36  000  and  the  momentum  thickness  Reynolds  number 


was  1780. 
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Table  3.1  Experimental  conditions  for  water  flow  with  no  injection. 


Um  =  0.574  m/s 

h  =  60.0xlCT3  m 

U0  *0.651  m/s 

v  =  0.957xl(T6m2/s 

Ut  =  27.6xl0“3  m/s 

Umh 

Re  =  =  36  000 

V 

5*  =  3.45xl(T3  m 

U0h 

Reg  =  °  =20  400 

2v 

0  =  2.61xl<T3  m 

U0h 

Ree  =  — =  1770 

The  shear  velocity  Ut  for  the  water  flow  with  no  injection  was  determined  using 
the  linear  variation  of  total  shear  stress.  Figure  3.1  shows  the  variation  of  calculated 
shear  velocity  with  distance  above  the  wall.  The  value  of  the  shear  velocity, 
determined  by  averaging  these  values,  was  0.0276  ±0.0010  m/s  at  95  per  cent 
confidence.  This  is  equivalent  to  a  strain  rate  of  798  s_1  at  the  wall  and  a  viscous  length 
v/Uj  of  0.0346  mm.  Therefore,  y+  =  867  at  the  channel  centerline. 

Since  the  flow  is  fully  developed,  static  pressure  drop  should  be  proportional  to 
the  wall  shear  stress  (if  one  assumes  that  this  high-aspect-ratio  channel  approximates 
flow  between  infinite  parallel  plates).  The  value  of  shear  velocity  calculated  from  the 
static  pressure  drop  is  about  ten  percent  higher  than  the  value  determined  from  the  total 
shear  stress  due  to  the  presence  of  secondary  flows  in  the  channel  comers. 


Shear  Velocity  (m/s) 
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3.2  Mean  velocities 

Figure  3.2  shows  the  mean  stream  wise  velocity  profile  normalized  with  shear 
velocity  and  kinematic  viscosity  (this  normalization  is  hereafter  indicated  by  the  + 
superscript).  The  two  solid  lines  in  the  figure  represent  the  near-wall  variation  in  U*. 

ir=y+'  (3.1) 

and  the  relationship  for  the  logarithmic  overlap  region 

U+  *  1  In  y+  +  B.  (3.2) 

K 

There  is  good  agreement  between  the  line  and  the  data  in  the  linear  region  and  there  is 
good  correspondence  with  the  logarithmic  law  using  ie=0.41  and  B=6.0  (Luchik  and 
Tiederman,  1987). 

33  Second  moments 

Figure  3.3  shows  the  root-mean-square  (RMS)  velocity  fluctuation  level  for  the 
stream  wise  and  normal  velocity  components  as  a  function  of  distance  from  the  wall. 
These  results  show  a  peak  in  u'  of  2.9  at  about  y+  =  15  and  a  value  of  about  one  at  the 
channel  centerline.  The  RMS  of  the  normal  velocity  shows  a  broad  peak  of  about  1.2 
centered  near  y+  =  100  and  drops  to  slightly  less  than  one  at  the  channel  centerline. 
These  values  for  u'  and  v'  agree  well  with  Kreplin  and  Eckelmann  (1979)  but  are  about 
ten  percent  higher  than  those  reported  by  Wei  (1987),  Hussain  and  Reynolds  (1978)  and 
Kim  et  al.  (1987). 

The  effect  of  gradient  broadening  on  the  measured  u'  levels  due  to  the  finite  size 
of  the  measuring  volume  is  estimated  to  introduce  about  six  percent  error  at  y+  =  5  and 
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less  than  one  percent  error  for  y+  >  10  (see  Karpuk  and  Tiederman,  1976).  The  results 
of  Luchik  and  Tiederman  (1986)  show  that  for  laser  velocimeter  measurements  the  RMS 
levels  are  not  affected  by  spatial  averaging  over  the  span  wise  length  of  the  probe  as  is 
die  case  with  hot-wire  and  hot-film  devices.  Therefore,  since  the  uncertainty  in  u,  is 
small,  these  results  are  believed  to  accurately  represent  the  turbulence  structure. 

The  near-wall  behavior  of  the  RMS  velocities  is  shown  in  Figure  3.4.  The 
behavior  of  u'  is  as  expected  near  the  wall;  the  RMS  of  the  normal  velocity,  however, 
approaches  a  constant  value  as  y  approches  zero.  The  results  of  Finnicum  and  Hanratty 
(1985)  and  Kim  et  aL  (1987)  indicate  that  v'  varies  as  y2  at  the  wall.  The  present  results 
do  not  appear  to  be  consistent  with  this  trend  and  therefore  the  accuracy  of  the  v' 
measurements  for  y+  <  10  is  questionable. 

The  RMS  velocities  when  normalized  by  the  local  mean  velocity  (turbulence 
intensity)  are  shown  in  Figure  3.5  and  the  near-wall  variation  is  shown  in  Figure  3.6. 
Near  the  channel  centerline,  both  quantities  approach  about  4  to  5  percent,  as  is 
expected.  The  streamwise  turbulence  intensity  is  approaching  a  value  near  0.4  at  the 
wall  which  is  slightly  higher  than  Kim  et  al.  (1987)  (0.36)and  considerably  higher  than 
Hanratty  et  al.  (1977)  (0.3).  The  normal  turbulence  intensity  behaves  as  expected, 
except  for  y+  <  10. 

For  a  fully  developed  channel  flow,  the  total  shear  stress  (mean  viscous  shear 
stress  minus  Reynolds  shear  stress)  varies  linearly  across  the  channel 

ttoui  =  -  puv  =  Tw*u(l  -  y/a)  (3.3) 

dy 

The  turbulent  shear  stress  uv  is  shown  in  Figure  3.7  normalized  with  Ut;  also  shown  is  a 
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line  representing  the  linear  variation  in  total  shear  across  the  channel.  Away  from  the 
wall,  where  the  mean  viscous  shear  is  small,  the  agreement  between  the  line  and  the 
data  is  excellent.  Equation  3.3  can  be  solved  for  -puv 

-puv  =  TwaU(l  -y/a)-p-~^ .  (3.4) 

dy 

Figure  3.8  shows  the  same  data  on  a  semi-log  coordinates,  which  allows  a  more 
detailed  examination  of  the  near-wall  region,  along  with  a  line  representing  Equation 
3.4  (the  values  determined  from  subtracting  the  the  mean  viscous  shear  stress  from  the 
total  shear  stress).  There  is  good  agreement  between  these  two  values,  although  there  is 
some  scatter  in  the  data  near  the  peak  in  -uv. 

The  correlation  coefficient  R^y  is  a  measure  of  how  strongly  the  u  and  v 
fluctuations  are  correlated,  i.e.  the  degree  to  which  the  velocity  fluctuations  contribute 
to  momentum  transport.  If  u  =  kv  where  k  is  a  constant,  the  data  is  perfectly  correlated 
and  every  velocity  fluctuation  contributes  to  the  net  momentum  transport  (Ruv  =  1);  if 
the  fluctuations  in  u  and  v  are  independent,  Rw  is  zero.  The  peak  value  of  Ruv,  shown 
in  Figure  3.9  as  a  function  of  y+,  is  about  0.35  which  is  lower  than  that  seen  in  the 
numerical  work  of  Moin  and  Kim  (1982)  ("0.4  at  Re  =  24  300)  and  Kim  et  al.  (1987) 
(0.45-0.50  at  Re  =  5600).  The  results  of  Wei  (1987)  decrease  from  0.45  to  0.35  as 
Reynolds  number  increases  from  30  000  to  80  000.  Assuming  the  trend  toward  lower 
peak  values  with  increasing  Reynolds  number  is  correct,  the  present  results  are 
compatible  with  the  numerical  results  although  they  disagree  with  those  of  Wei.  The 
occurrence  of  a  local  peak  in  R^  near  y+  =  12  which  has  been  seen  in  all  the  published 
channel  flow  computations  (Moin  and  Kim,  1982;  Kim  et  al.,  1987;  Moser  and  Moin, 
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1987)  does  not  appear  in  the  present  experimental  results. 


3.4  Reynolds  stress  production 


One  can  derive  transport  equations  for  the  elements  of  the  Reynolds  stress  tensor 
Wy  Although  these  are  not  "conservation"  equations  in  the  strict  thermodynamic 
sense,  the  terms  of  the  equations  represent  contributions  to  the  maintenance  of  a  ujuj  in 
a  unit  control  volume.  The  general  form  of  these  equations  (see  e.g.  Bradshaw,  1978)  is 
shown  in  Equation  3.5. 
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Choosing  different  values  of  i  and  j  yields  equations  for  the  six  independent  Reynolds 
stresses.  The  left-hand  side  of  this  equation  represents  the  time  rate  of  change  and 
advection  of  a  given  Reynolds  stress.  These  are  balanced  by  the  terms  on  the  right- 
hand  side:  Term  I  represents  production  of  Reynolds  stress  by  interaction  of  the 
turbulent  motion  with  the  mean  rate-of-strain  field.  Term  n  accounts  for  the 
generation,  destruction  or  redistribution  between  components  by  means  of  pressure 
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fluctuations  and  is  generally  referred  to  as  the  pressure-strain  correlation.  The 
pressure-strain  correlation  causes  the  normal  stresses  to  become  more  isotropic;  the 

i 

small  magnitude  stresses  are  increased  at  the  expense  of  the  larger  stresses.  Terms  HI 

and  IV  represent  transport  through  turbulent  fluctuations:  m  accounting  for  transport 

by  velocity  fluctuations  and  IV  for  pressure  fluctuations.  Viscous  transport  of  Reynolds  j 

stress  is  represented  by  term  V  and  term  VI  represents  generation  or  destruction  due  to 

viscous  stress  fluctuations. 

i 

The  behavior  of  the  terms  in  Equation  3.5  is  presented  in  detail  by  Man  sour  et  al. 

(1988)  and  Spalart  (1988)  for  a  fully  developed  channel  flow  and  a  boundary  layer, 
respectively.  These  results  were  obtained  by  direct  numerical  simulation  of  turbulent  j 

flow. 

Equation  3.5  can  be  specialized  to  the  case  of  stationary,  fully-developed,  two- 
dimensional  (on  average)  flow.  For  this  case,  the  left-hand  side  of  the  equation  is  zero 

and  the  right-hand  side  is  simplified.  For  the  streamwise  normal-stress  u2  the  result  is 
Equation  3.6. 
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The  dominant  source  term  is  the  production  term  (I)  and  it  is  balanced  primarily  by 
destruction  due  to  viscous  stress  fluctuations  (VI)  and  redistribution  to  the  other  normal 
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stresses  through  the  pressure-strain  correlation  (II).  Equation  3.7  shows  the  resulting 
equation  for  the  wall-normal  normal-stress  v2. 


m 


IV 


p! 

dv 

p 

dy 

w 

VI 


r 

dv 

2 

i 

■  N 

dv 

2 

i 

dv' 

2 

dx 

•  4 

T 

dy 

»  4 

dz 

■  m 

(3.7) 


For  this  component  of  the  Reynolds  stress  tensor,  there  is  no  type  I  production  term. 
The  major  source  term  is  the  pressure-strain  correlation  (II)  which  increases  v2  and 
decreases  u2.  This  source  is  balanced  by  the  viscous  destruction  term  (VI).  Equation 
3.8  is  for  the  Reynolds  shear  stress  uv. 
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As  in  Equation  3.6,  the  dominant  production  term  is  term  I.  The  results  of  Spalart 
(1988)  and  Mansour  et  al.  (1988)  show  that  in  the  shear  stress  equation,  the  pressure- 
strain  correlation  (II)  is  primarily  responsible  for  the  destruction  of  shear  stress  and  that 
destruction  through  viscous  stress  fluctuations  (term  VI)  is  very  small.  For  all 
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components  of  the  Reynolds  stress  tensor,  the  contributions  from  turbulent  diffusion 
(terms  III  and  IV)  and  viscous  diffusion  (term  V)  are  generally  small.  The  equation  for 

w2  is  similar  to  Equation  3.7  and  the  other  two  shear  stresses,  uw  and  vw  are  identically 
zero. 


Examination  of  Equations  3.6  and  3.8  shows  that  there  are  two  quantities 
responsible  for  the  production  of  Reynolds  stress  in  this  flow.  In  the  transport  equation 
for  the  streamwise  normal  stress  u2  (Equation  3.6),  the  production  term  has  the  form 


—  dU 
uv— —  • 
dy 

For  the  Reynolds  shear  stress  uv,  the  production  term  (from  Equation  3.8)  is 


(3.9) 
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Figure  3.10  shows  the  u2  production  profile  for  this  flow  normalized  with  shear  velocity 
and  kinematic  viscosity.  Also  shown  is  a  line  representing  this  quantity  calculated 
using  the  value  of  uv  from  Equation  3.4.  The  agreement  between  the  two  results  is 
excellent  and  they  compare  favorably  with  results  presented  in  the  literature.  Reynolds 
shear  stress  (uv)  production  is  shown  in  Figure  3.1 1.  This  quantity  exhibits  a  peak  near 
y+  =  10  which  is  significantly  higher  than  that  reported  by  Spalart  (1988),  Man  sour  et 
al.  (1988)  and  Nagano  and  Hishida  (1985).  Those  studies  both  showed  a  peak  of  about 
0.1-0.15  at  y+  =  15-20.  For  y+  >  10  the  behavior  is  reasonable  and  the  difference 


nearer  the  wall  is  attributable  to  the  error  in  the  measured  v'. 
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CHAPTER  4  -  THE  CONCENTRATION  FIELD  IN  INJECTED  FLOWS 

Concentration  measurements  were  made  while  injecting  either  pure  water  or  a 
700  ppm  aqueous  solution  of  SEP  ARAN  ap-273.  The  flow  in  the  channel  was 
maintained  at  the  conditions  described  in  the  previous  section  for  the  flow  with  no 
injection.  The  fluids  were  injected  at  both  the  top  and  bottom  channel  walls  at  a  flow 
rate  equal  to  the  mass  flow  rate  through  the  linear  portion  of  the  viscous  sublayer  of  the 
undisturbed  channel  flow  (1.00/min  through  each  injection  slot).  For  the  polymer 
solution,  this  resulted  in  25  percent  drag  reduction  in  the  region  from  50  to  150  mm 
downstream  of  the  injector  and  a  peak  drag  reduction  level  of  44  percent  in  the  region 
from  150  to  250  mm  based  on  the  measured  change  in  the  static  pressure  drop. 

For  the  flow  with  water  injection,  concentration  measurements  were  made  at 
locations  10,  25,  50,  and  100  mm  downstream  of  the  injector.  The  first  one  hundred 
elements  of  the  photodiode  array  were  sampled  resulting  in  concentration  profiles 
which  covered  6  mm  of  the  flow  starting  at  the  lower  channel  wall.  For  polymer 
injection,  concentration  measurements  were  made  at  locations  10,  25,  50,  100  and  200 
mm  downstream  of  the  injector  (x+  =  290,  730,  1450,  2910  and  5810  respectively). 
The  first  one  hundred  elements  of  the  photodiode  array  were  sampled  and  the 
measurements  were  made  over  7  mm  adjacent  to  the  channel  wall. 
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4.1  Time-resolved  concentration  profiles 

In  order  to  examine  the  type  of  instantaneous  behavior  in  the  concentration  field 
which  gives  rise  to  the  time  averaged  statistics,  time-resolved  measurements  of 
injectant  concentration  profiles  are  presented  in  Figures  4.1  through  4.6. 
Concentrations  are  normalized  with  the  injected  concentration  and  the  distance  from  the 
wall  is  normalized  with  the  shear  velocity  and  kinematic  viscosity  of  the  water  flow 
with  no  injection.  These  data  were  taken  at  x  =  25  mm  for  the  flow  with  dyed  water 
injection  and  for  all  five  measurement  stations  in  the  flow  with  polymer  injection. 
Figures  4.1  through  4.4  are  representative  of  the  entire  time  records  from  which  they 
were  taken.  At  the  two  downstream  measurement  stations  there  were  short  periods  of 
activity  interspersed  with  considerably  longer  intervals  which  were  relatively  quiescent. 
(Here,  activity  is  defined  as  movement  of  high  concentration  fluid  away  from  the  wall.) 
Figures  4.5  and  4.6  were  taken  from  active  periods  and  are  not  representative  of  the 
entire  time  sequence  from  which  they  were  taken. 

For  the  flow  with  water  injection,  Figure  4.1,  there  are  periods  of  activity  and 
quiescent  periods  interspersed  in  time  and  relatively  high  concentrations  extending 
instantaneously  to  y+  =  150  and  beyond.  In  Figure  4.2  time-resolved  concentration 
profiles  at  x  =  10  mm  for  polymer  injection  are  shown.  The  concentration  at  the  wall  is 
consistently  high,  but  in  the  region  from  y+  of  25  to  100  the  concentration  varies 
between  one  and  zero.  This  indicates  that  the  initial  mixing  of  the  polymer  solution 
occurs  through  stretching  and  distortion  of  the  injected  polymer  layer,  rather  than 
through  diffusive  transport  The  regions  of  high  concentration  fluid  extending  from  the 
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wall  appear  to  persist  for  the  same  duration  as  the  periods  of  activity  in  the  water  flow. 
This  may  mean  that  the  timescales  associated  with  the  initial  mixing  of  the  polymer 
solution  with  the  flow  are  similar  to  those  of  the  event  responsible  for  transport  in  the 
flow  without  polymer  injection. 

At  x  =  25  mm  and  50  mm  (Figures  4.3  and  4.4)  the  concentrations  are  reduced 
and  the  flow  is  relatively  quiescent  with  very  little  movement  of  the  high  concentration 
fluid  outward  from  the  wall.  One  would  expect  very  little  change  in  the  long  time 
average  statistics  in  this  region  and  also  very  small  values  of  vc. 

The  last  two  stations  exhibit  considerably  more  high  concentration  fluid  moving 
away  from  the  wall.  The  spatial  behavior  of  the  concentration  field  can  be  inferred 
from  its  temporal  behavior  by  assuming  the  concentration  field  is  "frozen"  and  moves 
past  the  probe  at  a  fixed  convection  velocity  (i.e.  the  Taylor  hypothesis).  The 
trajectory,  in  time,  of  the  high  concentration  regions  indicates  that  high  concentration 
fluid  lifts  away  from  the  wall  in  filaments  which  rise  near  their  downstream  end  and 
remain  near  the  wall  at  their  upstream  end.  At  the  last  location  (Figure  4.5)  the 
concentration  levels  are  generally  lower  than  at  x  =  100  mm  (Figure  4.6)  and  there  is 
slightly  more  activity. 

4.2  Probability  density  functions  for  concentration 

Plots  of  the  probability  density  of  measured  polymer  concentration  as  a  function 
of  normalized  concentration  and  distance  from  the  wall  for  the  various  streamwise 
locations  are  shown  in  Figures  4.7  through  4.11.  The  plots  span  the  region  from  the 
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wall  to  where  there  is  minimal  probability  of  non-zero  concentration.  One  apparent 
anomaly  is  present  in  Figure  4.7,  the  plot  for  x  =  10  mm:  immediately  adjacent  to  the 
wall  (y  £  0.25  mm,  y+  £  7.3)  the  probability  of  occurrence  of  concentrations  near  the 
injected  concentration  (C/Q=l)  is  reduced  compared  to  larger  y-locations.  It  is 
believed  that  this  is  due  to  the  slight  difference  in  index  of  refraction  between  the  water 
and  polymer  solution  and  the  high  polymer  concentration  present  at  the  wall.  This 
apparent  inaccuracy  in  the  measured  concentration  near  the  wall  would  introduce  less 
than  two  percent  error  in  the  concentration  values  determined  for  locations  farther  from 
the  wall.  This  effect  was  not  apparent  at  any  of  the  other  streamwise  locations. 

At  x  =  10  mm  (Figure  4.7),  the  concentration  at  the  wall  varies  from  about  0.35 
up  to  1.0.  At  about  y+=25,  the  presence  of  pure  water  is  detected  and  at  this  location  all 
concentrations  between  zero  and  1.0  were  measured.  This  behavior  will  result  in  large 
RMS  concentration  levels  at  this  location.  The  measured  concentrations  steadily 
decrease  with  distance  from  the  wall  until  virtually  no  polymer  is  present  at  about 
y+  =  80.  At  x  =  25  mm,  shown  in  Figure  4.8,  the  distribution  of  concentrations  near 
the  wall  has  broadened  to  include  concentrations  less  than  0.2.  Pure  water  is  still  not 
detected  for  y+  <  25  and  the  distributions  at  all  y-locations  except  those  immediately 
adjacent  to  the  wall  are  becoming  narrower  and  shifting  toward  lower  concentrations. 
The  probability  density  function  at  x  =  50  mm,  shown  in  Figure  4.9,  is  virtually 
unchanged  from  that  at  x  =  25  mm  and  is  further  evidence  that  little  transport  away 
from  the  wall  has  occurred  over  this  distance. 
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Figures  4.10  and  4.11  show  the  last  two  streamwise  stations.  At  these  locations, 
the  distributions  near  the  wall  become  less  broad  and  move  toward  zero  concentration 
while  presence  of  the  polymer  solution  is  detected  at  increasing  distances  from  the  wall. 

It  is  interesting  to  note  that  at  no  time  is  the  presence  of  pure  water  detected  at  the 
wall.  This  indicates  that  transport  in  this  region  is  dominated  by  molecular  diffusion 
and  small  scale  turbulent  mixing  and  that  fluid  from  the  outer  flow  rarely  if  ever 
penetrates  all  the  way  to  the  wall. 

4.3  Time-averaged  statistics 

Figure  4.12  shows  mean  concentration  profiles  as  a  function  of  distance 
downstream  of  the  injector  for  injection  of  dyed  water.  Again,  the  profiles  are 
normalized  with  the  injected  concentration,  Q,  and  are  plotted  versus  distance  from  the 
wall  normalized  with  the  kinematic  viscosity  and  shear  velocity  of  the  water  flow 
without  injection.  At  x  =  10  mm,  the  mean  concentration  at  the  wall  is  reduced  to  0.65. 
The  near-wall  concentration  is  reduced  by  about  one  half  at  each  succeeding 
streamwise  location  until  it  reaches  less  than  0.06  at  x  =  100  mm  (when  the  polymer  is 
mixed  uniformly  with  the  flow,  the  concentration  is  0.0018).  The  thickness  of  the 
high-concentration  wall  layer  increases  steadily  with  streamwise  distance  from  y+  =  70 
at  x  =  10  mm  to  y+ >  150  at  x=100mm.  Figure  4.13  shows  mean  concentration 
profiles  as  a  function  of  distance  downstream  of  the  injector  for  SEPARAN  AP-273 
injection.  The  near-wall  mean  concentration  is  0.75  at  x  =  10  mm  although  this  is 
probably  slightly  low  due  to  the  apparent  concentration  measurement  error  near  the 
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Figure  4.12  Mean  concentration  profiles  for  water  injection:  a)  x=10mm,  b) 
x  =  25  mm,  c)  x  =  50  mm,  d)  x  =  100  mm. 
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wall  at  this  streamwise  location.  Near- wall  concentration  then  levels  off  at  about  0.62 
from  25  to  50  mm  downstream  of  the  injector.  At  x  =  100  mm  the  near-wall 
concentration  is  reduced  to  0.5  and  at  200  mm  it  is  down  to  0.25.  For  x  £  100  mm  the 
thickness  of  the  high-concentration  layer  at  the  wall  is  about  75  wall  units  and  only  at 
x  =  200  mm  does  it  extend  beyond  y+=100. 

Although  the  water  mixes  more  readily  with  the  flow  than  the  polymer  solution, 
the  mean  concentration  profiles  at  x  *  10  mm  are  similar  enough  to  indicate  that  the 
initial  mixing  is  not  strikingly  different  for  the  two  cases.  At  x  =  50  mm  the  mean 
concentration  profile  for  polymer  injection  is  virtually  unchanged  from  that  of 
x  =  25  mm.  This  is  evidence  of  a  polymer-induced  supression  oi  turbulent  mass 
transport  in  this  streamwise  region.  At  x  =  100  mm,  the  suppression  of  the  mass 
transport  has  decreased  as  indicated  by  the  reduction  in  wall  concentration  from  that  of 
the  previous  location  and  at  the  last  measurement  station  the  mean  concentrations  are 
reduced  considerably  and  the  thickness  of  the  high-concentration  wall  layer  is  growing. 

Figure  4.14  shows  the  root-mean-square  (RMS)  concentration,  c\  as  a  function  of 
y+  for  the  injection  of  dyed  water.  At  x  =  10  mm  the  c'  profile  peaks  at  about  0.3Q  at 
y+=12  and  fluctuations  exist  out  to  y+=100,  significantly  beyond  the  edge  of  the  high- 
concentration  wall  layer  (based  on  the  mean  concentration  profile).  With  increasing 
streamwise  distance,  the  peak  in  the  RMS  concentration  profile  broadens  and  moves 
further  from  the  wall  and  non-zero  RMS  levels  occur  well  beyond  the  y-location  where 
the  mean  concentration  profile  approaches  zero.  RMS  concentration  profiles  for  polymer 
injection  are  presented  in  Figure  4.15.  At  x  =  10  mm,  the  RMS  concentration  level  is 
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Figure  4.15  Root-mean-square  concentration  profiles  for  polymer  injection:  a) 
x  =  10  mm,  b)  x  =  25  mm,  c)  x  =  50  mm,  d)  x  =  100  mm,  e)  x  =  200  mm; 
- ,  water  injection  results. 
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less  than  that  for  the  water  injection  case  across  the  entire  flow.  This  is  indicative  of  the 
decreased  mixing  when  polymer  solutions  are  injected.  At  the  next  streamwise 
location,  x  =  25  mm,  the  peak  fluctuation  level  has  decreased.  The  reduced  fluctuation 
level  near  the  wall  shows  that  the  concentration  in  the  high-concentration  layer  is  more 
uniform  than  at  the  upstream  station.  Here,  the  fluctuation  level  approaches  zero  at  the 
edge  of  the  high-concentration  layer  indicating  that  there  is  minimal  mixing  between 
the  wall  layer  and  outer  region.  At  x  =  50  mm  the  RMS  concentration  profile  is 
essentially  unchanged  but  at  x  =  100  mm  and  200  mm  non-zero  fluctuation  levels 
extending  beyond  the  high-concentration  wall  layer  are  evidence  of  increased  transport 
of  the  polymer  solution  from  the  wall  region.  This  increased  transport  results  in  the 
decrease  in  near-wall  concentration. 


CHAPTER  5  -  THE  STRUCTURE  OF  TURBULENCE  IN  INJECTED  FLOWS 


To  examine  the  effect  of  the  evolving  polymer  concentration  field  on  the 
structure  of  turbulence,  two-component  velocity  measurements  were  made  during 
polymer  injection  at  the  locations  25,  50,  100  and  200  mm  downstream  of  the  injector. 
To  determine  the  effect  of  injection  alone,  two-component  velocity  measurements  were 
made  at  x  =  25  mm  in  the  flow  with  water  injection.  Unless  otherwise  indicated,  the 
results  are  normalized  with  the  shear  velocity  of  the  water  flow  with  no  injection.  The 
distance  from  the  wall  is  normalized  with  this  shear  velocity  and  the  kinematic 
viscosity  of  the  channel  water.  These  normalizations  are  for  reference  only  and  are  not 
meant  to  imply  any  scaling  arguments.  Included,  also  for  reference  are  lines 
representing  the  appropriate  results  from  the  water  flow  presented  in  Chapter  3. 

5.1  Mean  velocities 

The  mean  velocity  profile  for  the  flow  with  water  injection  is  shown  in  Figure 
5.1.  There  is  a  slight  acceleration  near  the  wall,  and  a  slight  deceleration  for  y+  from  30 
to  about  100.  Outside  of  y+  of  200  the  injection  has  negligible  effect.  The  region 
where  the  flow  is  modified  corresponds  roughly  to  the  portion  of  the  flow  where  the 
mean  concentration  is  non-zero.  Figure  5.2  shows  the  mean  velocity  profiles  for  x- 
locations  between  25  and  200  mm  downstream  of  the  injectors  during  polymer 
injection.  In  contrast  to  the  results  for  water  injection,  the  outer  flow  (y+  >  50)  has 


accelerated  slightly  for  all  streamwise  locations  while  near  the  wall  the  flow  has 
decelerated.  This  is  the  opposite  of  what  occurred  for  water  injection.  The  region  of 
deceleration  corresponds  to  the  high-concentration  wall  layer  indicating  that  the 
deceleration  is  related  to  the  presence  of  the  polymer  solution  (due  to  both  an  increase 
in  viscosity  and  the  onset  of  drag  reduction).  The  reduction  in  velocity  near  the  wall, 
which  is  responsible  for  the  acceleration  of  the  outer  flow,  is  relatively  constant  for  the 
first  three  locations  and  is  slightly  smaller  at  the  200  mm  station. 

52  Second  moments 

Root-mean-square  velocity  fluctuations  normalized  with  the  shear  velocity  for  the 
flow  with  no  injection  are  shown  in  Figure  5.3.  It  is  clear  from  this  figure  that  while  u' 
is  virtually  unchanged,  the  RMS  of  the  normal  velocity  is  increased  across  the  high- 
injectant-concentration  portion  of  the  flow  compared  to  the  flow  with  no  injection.  The 
turbulent  shear  stress  for  the  water  injection  flow  is  shown  in  Figure  5.4.  Compared  to 
the  water  flow,  there  is  a  general  increase  in  -uv  in  the  near-wall,  high  concentration 
region.  These  results  show  that  the  major  effects  of  the  injection  process  on  the 
structure  of  turbulence  are  to  increase  the  turbulent  shear  stress  and  the  RMS  of  the 
normal  velocity  in  the  near-wall,  high-injectant-concentration  region. 

For  polymer  injection,  the  RMS  fluctuation  levels  for  the  streamwise  velocity  are 
shown  in  Figure  5.5.  There  is  a  15  percent  increase  in  peak  fluctuation  level  at  all  x- 
locations  and  the  peak  has  moved  outward  from  y+  =  15  to  about  30.  Near  the  wall,  the 
RMS  level  has  decreased  while  in  the  region  where  the  polymer  concentration  is  small 
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Figure  5.5  Root-mean- square  velocity  fluctuation  (u'/uT)  profiles  for  flow  with  polymer 
injection:  a)  x  =  25  mm,  b)  x  =  50  mm,  c)  x  =  100  mm,  d)  x  =  200  mm; 
- ,  no-injection  results. 
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(y+  £  50)  there  is  no  change  relative  to  the  water  flow.  Figure  5.6  shows  the  RMS 
fluctuation  level  of  the  normal  velocity.  At  x  =  25  mm  the  RMS  level  is  increased  near 
the  wall,  similar  to  the  behavior  in  the  flow  with  water  injection.  Near  the  wall,  v' 
decreases  with  streamwise  distance  except  for  y+  <  10  where  there  is  no  change.  (Due 
to  the  apparent  inaccuracy  of  the  measurement  of  v'  in  the  region  y+  <  10  it  is  likely 
that  this  does  not  accurately  reflect  the  turbulence  structure  for  y+  <  10.)  The  damping 
of  v'  extends  significantly  beyond  the  high-concentration  wall  layer  (based  on  the  mean 
concentration  profile)  for  x  =  100  mm  and  200  mm.  Similar  behavior  was  seen  by 
Walker  and  Tiederman  (1987).  This  previous  study,  which  examined  only  mean 
concentration  profiles,  did  not  show,  however, that  the  polymer  solution  instantaneously 
extends  beyond  y+  =  150  yielding  non-zero  RMS  concentration  levels  in  this  region. 

Reynolds  shear  stress  profiles  for  polymer  injection  are  shown  in  Figure  5.7.  At 
x  =  25  mm  the  maximum  in  -uv  is  about  15  percent  higher  than  the  water  flow.  This 
level  is  significantly  higher  than  for  the  flow  with  water  injection.  Elevated  -uv  levels 
persist  to  the  next  streamwise  station  x  =  50  mm.  In  a  manner  similar  to  that  for  v',  uv 
is  reduced  in  the  near-wall  region  and  this  effect  extends  beyond  the  high-concentration 
layer  at  the  two  downstream  locations. 

The  correlation  coefficient  Ruv  is  shown  in  Figure  5.8  for  the  flow  with  polymer 
injection.  At  x  =  100  mm  and  200  mm,  Ruv  is  significantly  decreased  in  the  near-wall 
region.  This  means  that  the  decreased  -uv  levels  which  occur  at  these  streamwise 
locations  are  not  a  result  of  a  general  decrease  in  the  fluctuation  level,  but  result  from  a 
smaller  contribution  from  the  fluctuations  to  net  momentum  transport. 
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S3  Reynolds  stress  production 

Figures  5.9  and  5.10  show  the  production  terms  which  appear  in  the  Reynolds 
stress  transport  equations.  For  the  normal  stress  u2,  shown  in  Figure  5.9,  production 
(term  I  in  Equation  3.6)  is  increased  at  x  =  25  mm  and  50  mm  due  to  the  increase  in  uv 
as  a  result  of  the  injection  process.  At  x  =  100  mm  the  production  has  significantly 
decreased  in  the  high-mean-concentration  wall-region  and  at  x  =  200  mm,  production 
has  been  virtually  eliminated.  Shear  stress  production  (term  I  in  Equation  3.8),  shown 
in  Figure  5.10,  shows  a  marked  increase  at  the  two  upstream  locations  due  to  the 
increase  in  v'  resulting  from  the  injection  process.  At  the  two  downstream  locations, 
production  is  essentially  the  same  as  in  the  water  flow. 

These  results  are  surprising  in  light  of  the  behavior  of  the  Reynolds  stresses 
themselves.  Although  normal-stress  production  is  virtually  eliminated  by  x  =  200  mm, 
u'  remains  elevated  and  essentially  constant  with  respect  to  distance  from  the  injector; 
meanwhile,  v'  decreases  with  stream  wise  distance.  An  increase  or  decrease  in  the  rate 
of  the  destruction  due  to  viscous  stress  fluctuations  (term  VI  in  Equations  3.5  through 
3.8)  would  affect  all  stresses  more  or  less  equally  since  over  the  range  of  length  scales 
which  contribute  to  term  VI  the  turbulence  is  nearly  isotropic.  As  was  discussed  in 
Chapter  3  (and  shown  in  Spalart,  1988  and  Man  sour  et  al.,  1988),  energy  is  transferred 
from  u2  to  v2  through  the  pressure-strain  correlations  (term  II).  An  increase  in  u2  and 
decrease  in  v2  indicates  that  this  transfer  has  been  inhibited  by  the  presence  of  the 
polymer  solution.  This  assumes,  of  course,  that  neither  the  Reynolds  stress  transport 
equations  nor  their  interpretation  is  altered  significantly  by  the  polymer. 


Figure  5.9  Reynolds  stress  production  (u2  component)  for  flow  with  polymer  injection: 

a)  x  =  25  mm,  b)  x  =  50  mm,  c)  x  =  100  mm,  d)  x  =  200  mm; - ,  no¬ 

injection  results. 
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The  turbulent  shear-stress  was  observed  to  decrease  significandy  with  stream  wise 
position  while  the  production  rate  for  this  stress  was  essentially  unchanged  from  that  of 
the  water  flow.  Spalart  (1988)  and  Man  sour  et  al.  (1988)  showed  that  viscous  stress 
fluctuations  play  only  a  minor  role  in  the  destruction  of  the  turbulent  shear-stress;  the 
dominant  quantity  is  the  pressure-strain  correlation  (term  II).  This  indicates  that  the 
mechanism  which  apparendy  inhibits  the  transfer  of  energy  between  the  normal- 
stresses  may  also  cause  increased  destruction  of  the  Reynolds  shear-stress. 
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CHAPTER  6  -  TURBULENT  MASS  TRANSPORT 

Turbulent  mass  transport  (vc)  profiles  were  measured  at  the  same  four  stream  wise 
locations  in  the  flow  with  polymer  injection  and  at  x  =  25  mm  in  the  flow  with  water 
injection.  Data  were  acquired  for  y-locations  extending  from  the  edge  of  the  high 
mean-concentration  layer  to  as  near  the  wall  as  possible. 

The  present  results  represent  a  unique  physical  situation  when  compared  to  scalar 
transport  results  presented  in  the  literature.  Most  scalar  transport  investigations  involve 
heat  transfer  with  either  constant  wall  heat-flux  or  a  constant  wall  temperature  and 
Prandtl  numbers  near  unity.  Due  to  the  similarity  in  boundary  conditions  for  the 
momentum  and  energy  equations,  the  structure  of  the  velocity  and  temperature  fields 
are  similar.  This  results  in  similarity  between  the  v0  and  uv  profiles  (see  e.g.  Antonia  et 
al.,  1977  and  Hishida  et  al.,  1986).  In  this  flow,  the  imposed  zero-flux  boundary 
condition  at  the  wall  (dC/dy  =  0)  and  the  high  Schmidt  number  destroy  this  similarity. 

6.1  Time-averaged  turbulent  mass  transport 

Figure  6.1  shows  the  vc  correlation  as  a  function  of  y+  for  the  flow  with  water 
injection.  Normalization  is  accomplished  with  the  injected  concentration  and  the  shear 
velocity  of  the  water  flow  with  no  injection.  The  vc  profile  does  not  appear  to  be 
similar  to  the  uv  profile  for  this  flow.  There  is  a  sharp  peak  in  the  profile  at  about 
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y+  =  40.  This  peak  corresponds  roughly  with  the  inflection  point  in  the  mean 
concentration  profile  but  does  not  correlate  with  any  major  feature  of  the  RMS  velocity 
or  concentration  profiles. 

The  vc  product  for  the  flow  with  polymer  injection  is  shown  in  Figure  6.2.  As 
expected,  vc  levels  are  near  zero  at  x  =  25  mm  and  50  mm.  Due  to  the  small  levels  for 
vc  the  uncertainty  is  high  (about  30  percent  of  the  measured  value);  however,  the  sign 
of  the  vc  product  is  believed  to  be  correct.  The  interpretation  of  these  negative  values 
is  uncertain.  At  the  two  downstream  locations,  the  turbulent  mass  transport  levels  off 
with  a  peak  at  about  0.01.  The  behavior  of  vc  is  in  good  agreement  with  the  foregoing 
discussion  of  the  mean  and  RMS  concentration  profiles.  The  vc  levels  during  polymer 
injection  are  significandy  less  than  those  measured  during  the  injection  of  dyed  water. 
This  is  not  unexpected;  since  the  polymer  requires  considerably  more  streamwise 
distance  to  mix  completely  with  the  flow,  the  mass  transport  normal  to  the  wall  must  be 
less. 

The  correlation  coefficient  for  mass  transport,  Rvc  =  vc/v'c',  for  the  flow  with 
water  injection  is  shown  in  Figure  6.3.  The  Rvc  profile  peaks  at  about  0.15  in  the  region 
from  y+  of  30  to  100.  This  indicates  that  the  normal-velocity  and  concentration 
fluctuations  are  less  strongly  correlated  than  uv.  Rvc  for  the  flow  with  polymer 
injection  is  shown  if  Figure  6.4.  (Note  that  the  scale  for  Figure  6.4  is  different  than 
Figure  6.3.)  At  the  two  upstream  locations,  Rvc  is  less  than  0.05  indicating  that  what 
fluctuations  exist  contribute  little  to  the  mass  transport  process.  Rvc  increases  with 
streamwise  distance  until  values  comparable  to  those  in  the  water  flow  occur  at 
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x  =  200  mm,  the  last  stream  wise  measurement  location. 

These  results  show  that  very  little  net  transport  away  from  the  wall  occurs  in  the 
polymer  flow  for  x  <  100  mm  as  was  inferred  from  the  mean  concentration  profiles,  the 
tune-resolved  concentration  profiles  and  the  probability  density  functions  for 
concentration.  The  occurence  of  comparable  values  of  the  correlation  coefficient  for 
x  =  200  mm  in  die  polymer  flow  and  x  =  25  mm  in  the  water  injection  case  means  that 
the  same  amount  of  mass  transport  occurs  in  each  flow,  given  the  magnitude  of  the 
velocity  and  concentration  fluctuations.  This  implies  similarity  in  the  ma^s  transport 
processes  in  these  two  cases. 

6 2  The  detailed  structure  of  turbulent  mass  transport 

In  order  to  better  understand  the  relationship  between  the  normal-velocity  and 
concentration  fluctuations  and  the  nature  of  the  vc  product,  a  more  detailed  analysis  of 
the  structure  of  the  normal-velocity  and  concentration  fluctuations  is  necessary.  Figure 
6.5  shows  scatter  diagrams  of  velocity  and  concentration  fluctuations  at  x  =  25  mm  and 
y+  =  88.3,  32.4  and  18.4  for  water  injection.  The  location  y+  =  32.4  corresponds  to 
where  vc  is  at  a  maximum  and  the  locations  on  either  side  of  this  maximum  are  where 
vc  was  equal  to  one-third  the  maximum  value.  At  all  three  locations,  the  distinct  cut¬ 
off  at  a  fixed  negative  concentration  fluctuation  value  corresponds  to  zero  concentration 
and  the  concentration  distributions  are  highly  skewed  as  was  previously  shown  in  the 
concentration  probability  density  functions.  At  low  concentrations,  all  possible  values 
for  the  velocity  fluctuations  are  present.  For  y+  =  32.4,  positive  concentration 
fluctuations  are  more  strongly  correlated  with  positive  velocity  fluctuations  resulting  in 


asymmetry  of  the  distribution  about  the  v  axis.  At  y+  =  88.3,  the  fluctuations  in 
velocity  are  comparable  to  those  at  32.4  while  the  concentration  fluctuations  are 
reduced  considerably  resulting  in  a  decrease  in  vc.  Positive  concentration  fluctuations 
are  more  strongly  correlated  with  positive  velocities  than  negative  velocities  here  also. 
The  near-wall  location  exhibits  concentration  fluctuations  that  are  comparable  to  those 
at  the  vc  peak  but  reduced  velocity  fluctuations  result  in  a  decrease  in  vc. 

To  determine  how  the  normal-velocity  fluctuations  present  at  a  point  contribute 
to  the  net  mass  transport,  the  joint  probability  density  function  for  velocity  and 
concentration,  f(v,c),  defined  such  that 

M  M 

j  j  f(v,c)dcdv=l  (6.1) 

was  calculated.  The  first  moment  of  the  joint  probability  density  function  is 

oo  m 

vc  =  J  J  [v  c  f(v,c)]  dc  dv  .  (6.2) 

Integrating  Equation  6.2  with  respect  to  concentration  only  and  normalizing  with  vc 
results  in 

oo 

1  =  J  F(v)  dv  (6.3) 

where 

oo 

F(V)  =  -L  J  [vcf(v,c)]dc.  (6.4) 

VC  -oo 

The  function  F(v)  shows  the  relative  contribution  of  normal-velocity  fluctuations  of  a 


given  magnitude  to  vc. 
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Figure  6.6  shows  F(v)  for  the  scatter  diagrams  in  Figure  6.5.  At  y+  =  88.3  the 
mass  transport  is  primarily  due  to  positive  velocity  fluctuations,  i.e.  high  concentration 
fluid  is  being  carried  away  from  the  wall.  At  the  vc  peak  there  is  a  substantial 
contribution  to  transport  by  negative  velocity  fluctuations  along  with  a  slightly  larger 
contribution  from  positive  velocities.  This  indicates  that  the  transport  consists  of  both 
movement  of  high  concentration  fluid  away  from  the  wall  and  movement  of  low 
concentration  fluid  toward  the  wall.  This  behavior  is  repeated  at  y+  =  18.4  although  the 
distribution  is  narrower  due  to  the  decrease  in  the  velocity  fluctuation  levels. 

Figure  6.7  shows  scatter  diagrams  of  velocity  and  concentration  fluctuations  at 
x  =  200  mm  and  y+  =  109.2,  27.1  and  13.1  for  polymer  injection.  Again,  the  middle 
location  corresponds  to  where  vc  is  a  maximum  and  the  locations  on  either  side  are 
where  vc  was  equal  to  one-third  the  maximum.  At  low  concentrations,  all  possible 
values  for  the  velocity  fluctuations  are  present.  For  y+  =27.1,  positive  concentration 
fluctuations  are  correlated  with  positive  velocity  fluctuations  resulting  in  asymmetry  of 
the  distribution  about  the  v  axis.  The  near-wall  location  again  exhibits  concentration 
fluctuations  that  are  comparable  to  those  at  the  vc  peak  but  the  velocity  fluctuations  are 
considerably  reduced.  This  results  in  the  reduced  value  for  vc.  At  y+  =  109.2,  the 
fluctuations  in  velocity  are  larger  than  those  at  the  vc  peak  while  the  concentration 
fluctuations  are  reduced  considerably  resulting  in  a  decrease  in  vc.  Here  also,  positive 
concentration  fluctuations  are  more  strongly  correlated  with  positive  velocities  than 
negative  velocities.  At  the  two  outer  locations,  there  is  a  stronger  correlation  between 
large  negative  velocity  fluctuations  and  near-zero  concentrations  (negative 
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concentration  fluctuations).  At  all  locations  the  distributions  appear  more  triangular  in 
shape  with  large  positive  concentration  fluctuations  corresponding  to  small  magnitude 

I 

velocity  fluctuations  and  large  velocity  fluctuations  occurring  only  when  the  polymer 
concentration  is  near  zero.  This  is  further  evidence  that  the  polymer  damps  the  normal 
velocity  fluctuations  and  demonstrates  a  direct  link  between  the  presence  of  the  I 

polymers  and  a  reduction  in  the  normal  velocity  fluctuation. 

Examination  of  F(v)  (Figure  6.8)  shows  behavior  very  similar  to  that  of  the  water 
injection  case.  One  difference  is  apparent:  the  transport  is  concentrated  at  smaller 
velocities.  This  is  due  to  the  generally  lower  v'  values  present  in  the  polymer  injection 
case. 

For  comparison  purposes,  Figure  6.9  shows  scatter  diagrams  for  x  =  25  mm  with 
polymer  injection  at  locations  comparable  to  those  shown  in  6.5  for  water  injection. 

The  values  of  vc  at  these  locations  are  all  small  and  less  than  zero.  At  the  outer 
location  the  distribution  is  similar  to  that  of  the  water  flow,  but  the  two  inner  locations 
exhibit  the  same  damping  of  the  normal  velocity  at  high  concentrations  exhibited  at 
x  =  200  mm.  These  distributions  have  values  of  v'  and  c'  comparable  to  those  seen  in 
the  near-wall  region  of  the  water  flow  but  the  distributions  are  nearly  symmetric  about 
the  v-axis  resulting  in  the  small  values  of  vc. 

Plots  of  F(v)  are  shown  in  Figure  6. 10  for  these  locations.  At  y+  =  32.4,  there  are 
contributions  from  both  positive  and  negative  velocity  fluctuations  that  result  in  net 
mass  transport  toward  the  wall.  A  contribution  from  positive  velocity  fluctuations  to 
mass  transport  toward  the  wall  is  displayed  at  y+  =  27.2  although  vc  is  very  small  here 
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and  there  is  significant  scatter  in  the  results.  The  location  at  y+  =  70.8  displays 
behavior  that  was  not  seen  in  the  previous  results:  The  contribution  from  positive 
velocity  fluctuations  is  of  opposite  sign  than  that  of  the  negative  velocity  fluctuations. 
This  indicates  that  at  this  location  (near  the  edge  of  the  high-mean-concentration  wall- 
region)  high  concentration  fluid  is  as  likely  to  be  moving  instantaneously  toward  the 
wall  as  away  from  the  wall.  This  results  in  near  zero  net  mass  transport 

In  the  flow  with  water  injection,  for  the  y-location  where  the  maximum  turbulent 
mass  transport  occurs  and  for  points  nearer  the  wall,  the  positive  and  negative  wall- 
normal-velocity  fluctuations  contribute  equally  to  mass  transport  indicating  that 
turbulent  transport  is  about  equally  divided  between  high  concentration  fluid  moving 
away  from  the  wall  and  low  concentration  fluid  approaching  the  wall.  Further  from  the 
wall,  mass  transport  is  primarily  due  to  positive  velocity  fluctuations  carrying  high- 
concentration  fluid  outward;  negative  velocities  carry  as  much  high-concentration  fluid 
as  low-concentration  fluid  inward  toward  the  wall  resulting  in  no  net  contribution  to 
mass  transport.  This  behavior  is  echoed  at  x  =  200  mm  for  the  polymer  flow  and  is 
evidence  that  the  transport  process  is  similar  for  these  two  flows.  At  x  =  25  mm  in  the 
flow  with  polymer  injection,  far  from  the  wall,  the  contribution  from  positive  velocities 
is  balanced  by  a  negative  contribution  from  negative  velocity  fluctuations.  Nearer  the 
wall,  positive  velocity  fluctuations  contribute  to  the  net  transport  but  negative 
fluctuations  as  a  group  do  not. 

The  previous  paragraph  concerns  net  transport.  In  the  scatter  diagrams  shown  in 
Figures  6.5,  6.7  and  6.9,  the  relative  symmetry  of  the  distributions  about  the  c-axis 
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accounts  for  the  low  values  of  vc  and  Rvc  in  these  flows.  This  is  probably  a  high 
Schmidt  number  effect  If  the  Schmidt  number  is  high  (the  mass  diffusivity  is  small 
relative  to  viscosity),  a  fluid  element  with  high  concentration  and  low  momentum 
which  moves  from  the  near  wall  region  to  the  low-concentration  outer  flow  retains  its 
high  concentration  even  after  it  has  come  to  equilibrium  with  the  momentum  field  in 
the  outer  flow.  When  this  fluid  element  is  carried  back  into  the  near-wall  region  (as  it 
inevitably  will  be)  it  returns  with  the  higher  momentum  of  the  outer  flow  but  still 
retains  the  high  concentration  it  had  when  it  left  the  near-wall  region.  This  scenario 
results  in  significantly  less  mass  transport  than  momentum  transport  ocurring  for  the 
same  turbulent  motion  and  explains  the  degree  of  symmetry  in  the  scatter  diagrams.  It 
also  explains  the  negative  contribution  to  mass  transport  from  large  negative  velocities 
at  locations  near  the  wall. 

63  Implications  for  modeling 

The  following  paragraphs  briefly  examine  the  nature  of  the  vc  product  and  its 
relationship  to  momentum  transport  for  the  purpose  of  modeling  this  type  of  flow. 
Comparisons  will  be  made  between  x  =  200  mm  for  the  flow  with  polymer  injection 
and  x  =  25  mm  for  the  flow  with  water  injection.  These  two  cases  exhibited  similar 
mean  conctration  profiles. 

Figure  6.11  shows  the  ratio  of  turbulent  mass  and  momentum  transport  (vc/uv) 
normalized  with  the  injected  concentration  and  the  shear  velocity  of  the  water  flow  with 
no  injection.  The  peak  in  this  ratio  for  the  polymer  injection  case  is  about  twice  that  for 
water  injection.  The  peak  for  water  injection  corresponds  to  the  peak  in  vc,  while  in  the 


polymer  injection  case  the  peak  is  nearer  the  wall  than  the  peak  in  vc.  These  results 
indicate  that  when  the  flow  is  modified  the  relationship  between  turbulent  mass  and 
momentum  transport  is  altered,  and  that  while  the  mass  and  momentum  transport  both 
decrease,  the  momentum  transport  is  reduced  more. 

Turbulent  diffusivities  are  defined  as 

vc 

dC/dy 

for  mass  transport,  and 

uv 

au/9y 

for  momentum  transport.  Profiles  of  turbulent  mass  and  momentum  diffusivities  are 
shown  in  Figure  6.12  normalized  with  the  kinematic  viscosity  of  the  water  flow.  The 
results  presented  are  for  x  =  25  mm  with  water  injection,  and  200  mm  for  the  flow  with 
polymer  injection.  These  data  show  that  the  mass-transport  diffusivities  are  very 
similar  for  the  two  flows.  The  momentum-transport  diffusivities,  however,  exhibit 
marked  differences.  The  values  for  the  polymer  flow  are  significantly  less  than  those 
for  the  water  flow.  This  is  further  evidence  of  the  difference  in  the  effect  of  the 
presence  of  the  polymer  solution  on  the  mass  and  momentum  transport  processes. 

Turbulent  Schmidt  numbers  are  shown  in  Figure  6.13  for  both  flows.  This 


quantity  is  defined  as: 


Sct  = 

Surprisingly,  the  flow  with  water  injection  has  a  higher  value  of  Set  than  the  flow  with 
polymer  injection.  The  high  values  of  Sct  exhibited  by  the  water  injection  case 
(Sct  £  1.5)  imply  that  the  accepted  values  of  0.7  -  1.0  are  inapproporiate  for  the  high 
molecular  Schmidt  number  regime.  The  low  values  exhibited  by  the  polymer  injection 
case  are  further  evidence  of  the  complex  relationship  between  mass  and  momentum 
transport  in  this  flow. 


123 


CHAPTER  7  -  SUMMARY  AND  CONCLUSIONS 

This  study  examined  a  fully  developed  turbulent  channel  flow  with  polymer 
injection  at  the  wall  and  had  two  complimentary  objectives.  The  first  objective  was  to 
gain  more  understanding  of  the  way  in  which  long-chain,  high  molecular  weight,  water 
soluble  polymers  cause  viscous  drag  reduction  in  turbulent  wall  flows.  The  second  was 
to  acquire  information  useful  in  the  modeling  of  the  transport  processes  in  flows  with 
polymer  injection  at  the  wall,  particularly  the  high  Schmidt  number  mixing  process 
which  occurs.  A  fluorescence  technique  was  developed  for  time-resolved 
measurements  of  species  concentration  profiles.  This  technique  was  combined  with 
laser  velocimetry  to  obtain  measurements  of  turbulent  mass  transport  vc  in  a  flow  with 
polymer  injection  and,  for  comparison  purposes,  a  flow  with  water  injection.  Two- 
component  laser  velocimeter  measurements  were  made  in  both  flows  to  examine  the 
effect  of  both  the  injection  process  and  the  evolving  polymer  concentration  field  on  the 
structure  of  turbulence. 

Two-component  velocity  measurements  were  used  to  examine  the  turbulence 
structure  in  the  basic  water  flow.  This  yielded  results  which  were  consistent  with  those 
in  the  literature  for  fully  developed,  two-dimensional  turbulent  channel  flow.  A  method 
was  presented  which  used  these  velocity  measurements  to  determine  the  wall  shear 
velocity  to  high  accuracy  (±3  percent). 
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Time-resolved  concentration  profiles  measured  during  polymer  injection  showed 
that  there  is  very  little  movement  of  high  concentration  fluid  away  from  the  wall  for 
x  <  100  mm.  This  was  indicative  of  the  absence  of  mass  transport  at  these  locations 
and  was  confirmed  by  the  behavior  of  the  mean  concentration  profiles  and  direct 
measurements  of  turbulent  mass  tansport.  For  x  >  100  mm,  there  were  periods  of 
activity  interspersed  with  somewhat  longer  quiescent  periods.  During  these  active 
periods,  long  (in  the  streamwise  direction)  filaments  of  high-concentration  fluid  moved 
away  from  the  wall;  this  is  the  dominant  transport  mechanism  in  the  near-wall  region. 
Compared  to  the  water  flow,  these  periods  of  activity  had  a  significantly  reduced 
frequency  in  the  polymer  flow. 

Examination  of  the  probability  density  functions  for  concentration  showed  that 
the  concentration  distributions  are  highly  skewed  toward  low  concentrations,  except 
very  near  the  wall.  They  also  revealed  that,  during  the  mixing  of  the  polymer  solution 
with  the  flow,  at  no  time  does  pure  water  penetrate  to  the  wall.  This  indicates  that  very 
near  the  wall,  the  transport  process  for  mass  (and  presumably  other  quantities)  is 
dominated  by  molecular  diffusion  and,  at  most,  small-scale  turbulent  mixing. 

A  detailed  examination  of  the  turbulence  structure  showed  that  in  the  flow  with 
polymer  injection,  the  RMS  streamwise  velocity  levels  are  increased  by  about  15  percent 
at  all  streamwise  locations.  Initially,  the  injection  process  causes  elevated  levels  of  the 
RMS  of  the  wall-normal  velocity  and  the  turbulent  shear  stress.  However,  damping 
action  by  the  polymer  solution  subsequently  reduced  these  quantities  significantly 
below  the  levels  present  in  the  water  flow  for  x  >  25  mm.  Scatter  diagrams  for  wall- 
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normal  velocity  and  concentration  showed  that,  instantaneously,  the  presence  of  high 
polymer  concentrations  resulted  in  small  velocity  fluctuations,  establishing  a  direct  link 
between  the  presence  of  the  polymer  solution  and  reduced  v'  values.  The  decreased 
-uv  levels  result  from  both  the  decrease  in  v'  and  a  decrease  in  the  degree  of  correlation 
between  the  u  and  v  velocity  fluctuations. 

An  examination  of  Reynolds  stress  production  showed  that  the  polymer  solution 
reduced  the  production  rate  for  u2  and  that  the  production  of  uv  was  not  changed.  The 
reduced  levels  of  v'  and  uv  along  with  the  increase  in  u'  caused  by  the  polymer 
indicates  that  the  polymer  may  alter  the  pressure- strain  correlations.  These  correlations 

transfer  energy  from  u2  to  v2  and  are  also  the  principal  source  of  uv  destruction. 
Assuming  that  the  Reynolds  stress  transport  equations  are  not  significantly  altered  for 
the  polymer  solution,  it  is  concluded  that  the  polymer  solution  inhibits  the  redistribution 
from  u'  to  v'  and  causes  increased  destruction  of  uv. 

In  general  vc  levels  were  lower  in  the  polymer  injection  flow  than  in  the  water 
flow.  This  results  in  the  polymer  requiring  a  longer  streamwise  distance  to  mix 
completely  into  the  flow  than  dyed  water.  The  correlation  coefficient  for  mass  transport 
was  lower  than  that  for  momentum  transport;  this  is  attributed  to  the  high  Schmidt 
number  associated  with  the  mass  transport  process.  The  low  molecular  mass  diffusivity 
(relative  to  viscosity)  results  in  high  concentration  fluid  mixing  from  the  wall  to  the 
outer  flow  and  back  to  the  wall  yielding  low  levels  of  net  turbulent  mass  transport 
relative  to  the  magnitude  of  the  velocity  and  concentration  fluctuations. 
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At  x=200  mm  in  the  polymer  flow,  vc  was  less  than  it  was  at  25  mm  with  water 
injection  but  the  mean  concentration  profiles  and  the  degree  of  correlation  were  similar 
in  both  cases.  This  indicates  that  the  reduced  transport  was  the  result  of  smaller 
velocity  and  concentration  fluctuations  and  not  differences  in  the  transport  process. 
The  contributions  to  mass  transport  from  the  wall-normal  velocity  fluctuations  were 
similar  in  these  two  situations.  A  brief  examination  of  modeling  issues  confirmed  these 
results:  Profiles  of  the  turbulent  diffusivities  for  mass  agreed  for  these  two  locations. 
However,  the  momentum  transport  was  significantly  reduced  in  the  polymer  flow 
causing  the  turbulent  Schmidt  number  profiles  to  be  different  in  these  flows.  For  these 
high  molecular  Schmidt  number  flows,  the  turbulent  Schmidt  numbers  deviated 
significantly  from  unity. 

From  the  results  of  this  study,  the  evolution  of  this  flow  from  a  Newtonian  state 
upstream  of  the  injector  to  a  drag-reduced  flow  can  be  described.  Immediately 
downstream  of  the  injector,  the  injected  polymer  layer  is  stretched  and  distorted  by  the 
water  flow  above  it.  The  larger  shear  viscosity  of  the  polymer  solution  causes  the  flow 
near  the  wall  to  be  slower  than  in  the  water  flow  with  no  injection.  This  blockage  of  the 
flow  causes  an  increase  in  the  velocity  gradient  at  the  edge  of  the  high  concentration 
layer.  This  along  with  the  disturbance  to  the  flow  from  the  injection  process  causes 
increased  Reynolds  stress  production  resulting  in  elevated  levels  of  v'  and  uv.  Since 
movement  of  high-concentration  fluid  from  near  the  wall  to  the  outer  flow  requires 
acceleration  in  the  direction  of  fluid  movement  (i.e.  an  elongational  strain)  the  high 
extensional  viscosity  of  the  polymer  solution  prevents  it  from  being  carried  into  the 
outer  flow.  (This  is  confirmed  by  the  strong  correlation  between  high  polymer 
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concentrations  and  small  wall-normal  velocities.) 

As  the  flow  proceeds  downstream,  the  concentrations  near  the  wall  are  reduced 
and  become  more  homogeneous  as  a  result  of  small  scale  mixing  within  the  near-wall 
layer.  There  is  still  no  transport  of  high  concentration  fluid  from  the  wall  but  the  effects 
of  the  injection  process  begin  to  diminish  as  v'  and  uv  levels  return  to  those  of  the  water 
flow  with  no  injection.  This  shows  that  above  the  high-concentration  layer  the  flow  is 
still  Newtonian  in  nature. 

As  the  concentration  in  the  near-wall  region  is  reduced,  the  extensional  viscosity 
is  reduced  and  high-concentration  fluid  begins  to  move  away  from  the  wall  in  the  form 
of  long  filaments.  Due  to  the  low  molecular  mass  diffusivity  of  the  polymer  solution, 
this  high  concentration  fluid  is  mixed  back  into  the  near-wall  region  resulting  in  small 
net  mass  transport  relative  to  the  size  of  the  velocity  and  concentration  fluctuations. 
The  high  extensional  viscosity  of  the  polymer  solution  causes  these  polymer  filaments 
to  resist  acceleration  and  inhibits  the  transfer  of  energy  from  u'  to  v'.  This  resistance  to 
acceleration  increases  the  local  instantaneous  shear-strains  which  causes  increased 
destruction  of  uv  through  the  pressure-shear-strain  correlation.  It  also  allows  this  fluid 
to  return  to  the  near-wall  region  with  a  smaller  than  expected  increase  in  momentum 
resulting  in  decreased  net  momentum  transport.  This  decrease  in  turbulent  transport 
results  in  reduced  viscous  drag. 
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Appendix  A:  Calculation  of  Rotated  Statistics 

The  ability  to  calculate  second  statistical  moments  in  a  coordinate  system  which 
is  rotated  relative  to  the  coordinate  system  in  which  the  measurements  were  obtained 
simplifies  processing  of  two-component  velocity  measurements.  Statistics  can  be 
calculated  on  the  measured  velocity  components  and  then  these  results  can  be  used  to 
calculate  the  desired  statistics.  This  eliminates  the  computational  expense  of  resolving 
the  individual  velocity  measurements  into  the  desired  components  prior  to  calculating 
the  statistics.  It  also  allows  one  to  change  the  coordinate  system  and  determine  the 
appropriate  statistics  without  re-analyzing  the  original  velocity  data. 

We  desire  to  calculate  statistics  in  the  x-y  coordinate  system  (U  and  V  velocity 
components)  using  statistics  for  Ui  and  U2  which  were  measured  in  the  xrx2 
coordinate  system  which  is  rotated  a  some  angle  0.  (See  Figure  A.l.)  The  relationship 
between  U  and  V  and  Ui  and  U2  is  expressed  as 

U  =  U2COS8  +  U 1  sin0  (A.l) 

V  =  U2sin0  -  Uicos0.  (A.2) 

Time-averaging  Equations  A.  1  and  A.2  yields 

U  =  U2COS0  +  Uisin0  (A. 3) 

V  =  U2sin0  -  Uicos0.  (A.4) 

Subtracting  the  equations  for  the  mean  velocities  from  those  for  the  instantaneous 
velocities  results  in  equations  for  the  fluctuating  velocities  u  and  v  in  terms  of  Ui  and 
U2.  These  equations  are 


1  Coordinate  system  for  rotation  of  velocity  statistics 
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u  =  U2COS0  +  Uj  sin0 

(A.5) 

v  =  U2sin0  -  U}cos0. 

(A.6) 

Squaring  Equation  A.5  and  time-averaging  results  in  an  equation  for  u2: 

u2  =  ui2sin20  +  U22cos20  +  2uiu2cos0sin0  (A.7) 

Performing  similar  operations  on  Equation  A .6  yields  an  equation  for  v2: 

v2  =  ui2cos20  +  U22sin20  -  2ui"u2cos0sin0  (A.8) 


Multiplying  Equations  A.5  and  A.6  and  time-averaging  yields  the  following  equation 
for  uv. 

UV  =  [u22  -  U! 2  j  sin0cos0  +  ujuj  £  sin20  -  cos20j  (A.7) 

These  equations  were  used  to  calculate  statistical  moments  in  the  desired  coordinate 
system  for  the  results  presented. 


Appendix  B:  Uncertainty  Analysis 


When  examining  profiles  of  turbulence  statistics  and  the  evolution  of  such 
profiles  in  the  direction  of  flow,  there  are  several  distinct  sources  of  uncertainty  to  be 
considered.  Absolute  errors  in  location,  as  well  as  errors  in  the  position  of  one  point 
relative  to  another  are  important.  In  the  estimation  of  turbulence  statistics  one  must 
consider  both  fixed  and  random  errors  in  the  individual  realizations  used  to  calculate  a 
statistic  and  the  inherent  uncertainty  in  estimating  a  statistic  with  a  finite-size  ensemble. 
This  section  is  an  attempt  to  evaluate  the  impact  of  these  uncertainties  on  the  results  to 
be  presented. 

Two  separate  traversing  mechanisms  were  used  to  move  the  laser  velocimeter 
measurement  volume  in  space.  The  vertical  (y-direction)  motion  was  controlled  by  a 
micrometer  having  a  least  count  of  0.013  mm  the  streamwise  (x-direction)  location  was 
determined  using  a  scale  with  a  least  count  of  0.5  mm.  This  results  in  probable  errors  in 
the  relative  locations  of  two  points  in  space  of  0.38  wall  units  in  the  y-direction  and  14 
wall  units  in  the  streamwise  direction.  The  absolute  location  of  the  wall  was 
determined  using  the  laser  velocimeter  measurement  volume  resulting  in  a  estimated 
uncertainty  of  0.9  wall  units  (one-half  the  y-extent  of  the  measurement  volume) 

Uncertainty  in  statistical  quantities  can  be  separated  into  two  distinct  quantities: 
the  probable  error  in  the  individual  realizations  used  to  calculate  the  statistics,  and  the 
uncertainties  associated  with  estimating  a  statistic  using  a  finite-size  ensemble.  After 
considering  each  of  these  in  turn,  this  section  closes  with  the  discussion  of  the 
combined  effect  of  both  types  of  uncertainty. 


The  results  of  this  study  are  derived  primarily  from  two-component  laser 
velocimeter  data,  one-component  laser  velocimeter  data  and  concentration  data 
acquired  using  the  laser-induced-fluorescence  technique.  For  the  laser  velocimeter 
data,  the  primary  source  of  error  is  the  uncertainty  in  the  determination  of  the  beam¬ 
crossing  angle  which  affects  the  calculation  of  the  fringe  spacing  This  results  in  an 
uncertainty  of  0.4  percent  for  the  measured  velocity  component.  (Since  the  statistics 
are  calculated  for  a  rotated  coordinate  system  where  V  =  0,  the  uncertainty  for  U  and 
for  V  is  about  0.004U.)  The  dominant  uncertainty  in  the  concentration  measurements 
results  from  the  pH  dependence  of  the  fluorescence  of  the  dye  used  and  the  uncertainty 
in  matching  the  pH  of  the  channel  water  to  that  of  the  injected  polymer  solution  (0.5 
percent  of  the  measured  concentration). 

These  errors  can  be  propagated  to  the  statistics  calculated  from  the  acquired  data. 
Experimental  uncertainties  (e’s)  for  the  mean  velocities  (in  the  x-y  coordinate  system) 
are 

eu,ev=  0.004  U;  (B.l) 

experimental  uncertainty  for  the  mean  concentration  is 

eg  =  0.005  C.  (B.2) 

For  the  RMS  quantities,  the  experimental  uncertainties  are 
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ec'  =  0.005  c'  (B.5) 

For  the  remaining  second  moments,  uv  and  vc: 

6(jv  =  0.004  £  +  uv  j  (B.6) 

Eve  ~  0.006  vc  .  (B.7) 

The  relatively  complex  form  of  the  uncertainties  for  uv  and  v'  results  from  measuring 
the  velocities  at  ±45°  to  the  desired  components. 

Appropriate  normalization  of  these  uncertainties  yields  : 


£u  £v 

— ,  —  =  0.004 , 

(B.8) 

u  u 

=  0.005, 

(B.9) 

C 

=  0.004 , 
u 

(B.10) 

^  =  0.004-"*-, 

V  v  2 

(B- 11) 

Ec' 

-V  =  0.005  , 
c 

(B.  12) 

=  0.004  f  1  +  ^-1  , 

(B.13) 

uv  [  UV  J 


and 
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—  =  0.006.  (B.14) 

vc 

The  estimation  of  any  statistic  from  a  finite-size  ensemble  of  individual 
realizations  has  an  inherent  uncertainty  which  is  defined  by  the  type  of  statistic  and  the 
size  of  the  ensemble  used  in  its  calculation.  The  standard  deviation  of  the  estimated 
mean  is 


1 

Ou  =  —f=  °u 
VN 


(B.15) 


If  one  assumes  that  the  distribution  for  a  given  quantity  is  Gaussian  or  nearly  so,  the 
uncertainty  at  the  95  percent  confidence  level  resulting  from  the  estimation  of  the 
statistic  from  a  finite  size  ensemble  (5)  is 


8p  _  2  u* 

u  €  u’ 


(B.16) 


This  indicates,  as  was  shown  by  Yanta  and  Smith  (1973),  that  the  uncertainty  in  the 
estimated  mean  depends  on  both  the  ensemble  size  and  the  RMS  level. 


Similar  uncertainties  can  be  calculated  for  the  second  moments.  For  large  sample 
sizes,  the  standard  deviation  of  the  standard  deviation  is  given  by 


<V  = 


>6n 


ay . 


(B.17) 


When  estimating  the  uncertainty  for  products  such  as  uv,  if  we  assume  the  two 
quantities  are  independent  then 
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0>uv  — 


CT|jOv  • 


The  uncertainties  for  u'  and  uv  are  therefore  given  by 

Su'  _  2 

u'  >ST  ’ 


and 


(B.18) 


(B.19) 
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uv 


uV 

uv 


_2 
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(B.20) 


where  R  is  the  correlation  coefficient  for  u  and  v.  This  shows  that  the  uncertainty  in  the 
estimate  of  an  RMS  quantity  depends  solely  on  the  ensemble  size  (Yanta  and  Smith, 
1973),  but  the  uncertainty  in  quantities  like  uv  depend  upon  both  the  ensemble  size  and 
the  degree  of  correlation. 


Figure  B.l  shows  the  variation  of  this  type  of  uncertainty  with  ensemble  size  and 
the  relevant  parameters  for  mean,  RMS  and  uv-type  statistics. 

The  total  uncertainty  in  a  given  statistic  is  then  the  root-sum-square  of  the 
uncertainties  arising  from  both  sources 

V  =  V52+e2  .  (B.21) 


For  the  mean  quantities: 
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For  the  second  moments: 
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These  equations  were  used  to  estimate  uncertainty  at  the  95  percent  confidence 
level  for  the  experimental  data.  For  mean  and  root-mean-square  statistics,  the 
uncertainty  estimate  is  generally  comparable  to  the  size  of  the  plotting  symbol  used. 
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For  uv  and  vc,  the  uncertainty  estimates  are  indicated  by  error  bands  on  a  few 
representative  points. 
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